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ABSTRACT 
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The  object  of  tbo  present  work  is  the  investigation  into  the 
probleas  of  stress  and  deforaation  of  linear  viscoelastic  solids 
under  various  environmental  conditions,  such  as  gravity,  temperature 
and  inertia  forces. 

J 

introduction,  which  deals  with  the  state  of  the 
art  fS-date,  the  body  of  the  report  is  divided  in  three  chapters. 

•in  the  first  chapter,  the  thermal  stresses  in  viscoelastic  solids 
vith  teaperature  dependent  properties  are  investigated,  teaperature 
dependence  being  limited  to  the  the  rmorheologl cally  simple  type.  /^Here 
approximate  methods  are  developed.  Finally,  the  infinite  cylinder 
under  a  transient  radial  teaperature  field  is  studied  and  two  analytic 
solutions  are  given 


^Tbe 


second  chapter  is  occupied  vith  the  evaluation  of  the  effects 
of  the  acceleration  terms  on  the  stress  distribution  in  incompressible 
viscoelastic  solids,  again  under  conditions  of  non-unifora  transient 
teaperature.vAnalytlc  solutions  are  given  to  the  probleas  of  the  hollow 
sphere  and  cylinder.  J 

CS»  third  chapter  is  devoted  to  gravitational  effects'.  Hie  problem 
of  the  horizontal  sluap  of  a  viscoelastic  hollow  cylinder  contained  in 
a  thin  elastic  shell  and  resting  on  a  rigid  horizontal  plane  is  solved 
by  the  application  of  the  Finite  Fourier  Transfora  and  the  Laplace 

Transform.  ____  _ 

*rS"  appendix  I  the  problea  of  an  infinite  slab  and  a  solid  sphere 


are  solved  nuatrlcally  for  the  purpose  of  illustration  and  compared 
with  exact  solutions.  n»  aethods  of  the  first  chapter  are  also  used 


for  numerical  solution  of  the  problem  of  cylinder  in  connection  with 
dig! tel  computer  progress.  In  Appendix  II  error  estimates  connected 
with  one  of  the  analytic  methods  axe  derived.  In  Appendix  III  the 
convergence  of  the  first  analytic  method  is  established. 
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GENERAL  IHTR0DUCT10M 

The  method  of  analysis  of  the  quasi- static  boundary  value  problem 
for  homogeneous  isotropic  linear  viscoelastic  materials  under  surface 
tractions,  body  forces  and  isothermal  temperature  conditions  is  now 
veil  established. 

The  formal  solution  of  the  above  problem  within  the  scope  of  small 
deformations  and  time  independent  boundaries  can  be  effected  by  the 
application  of  Laplace  transform  thereby  eliminating  the  time  depen¬ 
dence  and  thus  reducing  the  viscoelastic  problem  to  an  "associate" 
elastic  one. 

The  introduction  of  a  temperature  field  however,  be  it  transient 
or  steady,  introduces  new  difficulties  that  are  not  easily  surmountable- 

The  temperature  effects  are  two- fold,  i.e.  (a)  Thermal  strains 
are  set  up,  and  (b)  the  mechanical  properties  being  extremely  sensitive 
to  temperature  variations  -  orders  of  10°  C  are  significant  -  ensure 
that  a  non- homogeneous  transient  temperature  entails  a  non-homo geneoua 
material  with  mechanical  properties  as  functions  of  the  space  variables 
and  time. 

The  mode  of  variation  of  the  mechanical  properties  with  temperature 
is  in  itself  a  problem  which  has  only  been  partially  solved. 

In  the  present  paper,  these  difficulties  will  be  formally  investi¬ 
gated  and  some  methods  for  over costing  them  will  be  given. 

Review  of  the  Isothermal  ^uaslstatic  Problem 

Given  a  viscoelastic  body  B  with  surface  3  under  surface  tractions 

body  forces  f^  and  surface  displacements  u^( s)  which  may  possibly 
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depend  on  time ,  the  object  is  to  determine  the  stress  and  displacement 
distribution  within  the  body. 

The  viscoelastic  constitutive  relations  in  Cartesian  coordinates 

* 

can  always  be  put  in  the  fora 
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G^  is  the  relaxation  modulus  in  shear  and  Gg  is  the  relaxation  modulus 
in  dilatation. 

The  complete  solution  is  obtained  if  the  following  relations  are 
satisfied. 

Equilibrium  condition 


1 


Strain  displacement  relation 


(* **0 


Ul 


j 


(5) 


*  Repeated  indices  denote  suasntion. 

**  A  coma  followed  by  an  index  J  denote  differentiation  with  respect 
to  the  co-ordinate. 
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and  the  boundary  conditions: 


00 

1A-  ^  At . 


15,  > 

where  ^  t  are  prescribed  on  the  part  3^  of  the  boundary 


(6) 


<0 

T. 


(7) 


where  T  ^2)  are  prescribed  on  the  remaining  part  Sg  of  the  boundary  3- 
The  usual  approach  is  to  remove  time  dependence  by  applying 
Laplace  Transform. 

Then  Eq.  (l)  and  (2)  become 


MvtO -- 


(0) 

(9) 


for  an  initially  unstressed  and  unstrained  state. 


^i>i 

-  ii  =° 

(10) 

=  i  i  “<,i 

(11) 
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IX. t-  on  3^ 

(12) 
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-  (K) 

•t  =  x.; 


on  Sg 


(13) 


Eq.  (8)  to  (13)  clearly  correspond  to  an  elastic  problem  where 
the  dependent  variables  and  elastic  constants  are  functions  of  the 
Laplace  Transform  parameter  p. 

The  elastic  equivalents  of  (8)  and  (9)  are: 

S;;  =  2  S 
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=  5  f"T  fc 


If  a  solution  to  the  associate  elastic  problem  can  be  found  then 
the  viscoelastic  solution  is  derived  by  applying  the  inverse  Laplace 
Transform  to  the  elastic  solution  after  substitution  of  G^p  for  2G  and 
GgP  for  3K. 

The  method  of  solution  of  the  isothermal  problem  is  therefore 
clearly  defined. 


Influence  of  Temperature  on  the  Mechanical  Properties 

It  has  been  established  by  experiments  that  increase  in  temper¬ 
ature  accentuates  creep  and  relaxation  rates,  that  is,  strains  Increase 
and  stresses  decay  faster  with  time. 

The  direct  experimental  determination  of  the  dependence  of  these 
rates  on  the  current  temperature  is  an  exceedingly  difficult  task,  as 
one  can  appreciate. 
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However,  a  wide  class  of  morphous  polymers  obeying  the  linear 
viscoelastic  law  exhibit  a  simple  property  with  change  of  temperature. 
This  property  has  been  utilized  to  derive  the  dependence  of  the 
mechanical  properties  on  temperature  for  such  a  class.  This  property 
is  illustrated  in  Fig.  1. 


In  this  figure  tbe  relaxation  modulus  is  plotted  against  logt  for 
different  constant  uniform  temperatures. 

It  is  seen  that  the  effect  of  a  temperature  increase,  is  a  shift 
of  the  whole  relaxation  curve  to  the  left,  this  shift  being  a  function 
of  temperature. 

Let 


Jloo^ "t  -  /IL 


(lk) 


i 


M, 


~  -e 


(15) 


then  is  obtained  from  G(t)  by  the  transformation 


the  star  Is  used  to  indicate  the  change  in  functional  dependence. 
In  view  of  Fig.  1 


*  * 

G  (*)  =  Gj*<- 

T,  To 


(17) 


*  -ffn)  I 

G  Tlf*>-GT.f+«  j  <l8> 

|>(T„)  f(T)  >  O  T>  Te 

We  now  put 


e  =  afT,)  (19) 


Then 


(20) 


and  generally 
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We  nov  define  a  "reduced  tiae"  a  such  that 

f  =  a(T)  t  (22) 

Then  finally 

Q  (t)  =  G  (J)  (23) 

T  To 

Thus  knowing  G  at  some  uniform  reference  temperature  T0,  G  at  any 
other  uniform  temperature  can  be  found. 

From  (23) 


cL&t 

dt 


cL  6Ta((§)  _  d&r.  a(T)  (24) 

dg  dt  dt 


Eq.  (24)  gives  a  relation  between  the  relaxation  rates  at  temperature 
T  as  compared  to  To. 

For  Instance  at  t  =  0,  §  *  0 


d&T  I  _  clGr.  a  (T) 

cL“t  I  t  -o  d  ^  t*o 


(25) 


Note  that  |  ^ 
because  of  (17). 


for  all  temperature  histories  since 


ffr) 

cKt)s  e  ^  o 
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From  (22)  and  (23) 

Gr  CO  =  G_  cO  (26) 

lo 

Eq.  (26)  implies  that  the  elastic  (initial)  response  of  the 
material  is  unaffected  by  temperature  changes.  Also  in  view  of  (26) 
and  (23) 

(27) 


The  constitutive  relations  at  temperature  T  take  the  form: 


/ 

where  F  -  (X  ,  and  G,  is  referred  to  th  •  reference  temperature 

t 

To.  Similarly, 


3  fexk' 

9* 


‘(29) 


assuming  a  stress  free  field  in  the  presence  of  the  uniform  temperature 
T.  Experimental  observations  as  well  as  thermodynamic  considerations 
indicate  that  the  same  shift  factor  applies  both  to  and  Gg. 
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Steady  Non- uniform  Temperature 

Relations  (26)  and  (29)  can  be  ionediately  generalised  to  steady 
non- unifora  temperatures ,  where  now  both  a^  and,  hence  £  are 
functions  of  the  space  variables 

In  an  explicit  form  (28)  and  (29)  become: 

t 

S-  --  [  G  (30) 

o 

and 

t 

=  f  (31) 

"  o 

where 

T 

0  =  —  \  cJ  (j)  cLT  (32) 

o (0  J  — 

•o 

It  is  noteworthy  that  L.T.  is  still  applicable  to  (30)  and  ( 31) 

i.e. 


(33) 


* 


Henceforth  L.T.  will  signify  Laplace  Transform. 
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0.M  1 
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(3M 


However,  in  the  resulting  associate  elastic  problem  the  associate 
"elastic  constants"  are  nw  functions  of  the  apace  variables.  Closed 
fora  solution  to  such  a  problea  is  rather  unlikely,  except  for  the 
simplest  geometries,  and  resort  to  numerical  procedures,  is  inevitable. 


Temperature  as  Fun  tlon  of  Space  Variables  and  Time 

The  constitutive  relations  in  the  presence  of  transient  temper- 
ature  fields,  have  been  formulated  as  a  generalization  of  Eq  (28)  and 
(29) ,  and  take  the  following  fora 


where  now 


I  r  ^  &  lT g  (*,,«■)  (37) 

O 


*  See  Ref  [3] 
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Relations  (35)  and  (36)  though  still  relatively  siaple  sore  no 
longer  in  a  convolution  form  and  L.T.  is  inapplicable,  with  respect  to 
t. 

A  transformation  has  been  suggested  that,  apparently,  overcomes 
this  difficulty. 

Let 


f  -  f  Ok  .  ^  (39) 

(4o) 

It  must  be  noted  that  since  a<p  is  non- negative  ^(*<,0  is  a 
monotomlcally  increasing  function  of  time  for  any 

Hence  relation  (37)  may  be  inverted  in  the  fora: 


i 


(41) 


Then 


(42) 


and  similarly  for  other  functions. 

In  view  of  (39) »  (40),  (4l)  and  (42),  relations  (35)  and  (36) 


K 


l 

--  C  G(M') 


9€;; 


(43) 


take  the  fora: 


12 


Eq.  (43)  and  (W*)  are  now  convolution  integrals  in  the  ^ 
variable,  and  hence  L-T.  is  again  applicable  in  this  particular 
variable . 

This  transformation,  however,  modifies  the  field  equations  as  veil. 
For  the  sake  of  argument  let  a  function 


Then 


f  =  -f  (%  X  .  O 


ii  - 


2k 


(*5) 


(W) 


In  terms  of  the  new  variables 


>[  &  +  7>] 


( i  ■  sunned) 


or, 


•}X« 


t,x;( 


d- 


A 


(47) 


(U8) 


The  equilibrium  relation 
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*Li  »r  +ft 


-  o 


(49) 


becomes : 


4 

e, 


a  a.-, 

®T  *\ 


f;  =° 


(50) 


Note  that 


3f.  I  =■  F(x;,0  =  Hv  .  Hx  f; 


oi; 1 1 


(51) 


Since  L-T.  of  an  ordinary  product  of  two  functions  is  not  always  defined 
L.T.  is  inapplicable  to  Eq.  (50). 

Obviously  the  above  transformation  simplifies  the  constitive 
relations  but  complicates  the  field  equations. 

This  transformation,  however,  is  useful  in  the  case  where  T  *  T(t) 
only.  Then, 

f  (52) 

and  (50)  reduces  to 


A 


=  O 


(53) 


Now  L.T.  of  the  field  equations  and  constitutive  relation  can 
be  taken  with  respect  to  £  and  the  viscoelastic  problem  reduces  to 
an  associate  elastic  problem  as  in  the  isothermal  case. 

We  summarize  as  follows: 

(a)  Isothermal  case:  L.T.  is  applicable  with  respect  to  t  and  the 
viscoelastic  problem  reduces  to  an  associate  elastic  problem. 

(b)  Temperature  space  dependent  but  steady:  L.T.  applicable  with 
respect  to  t.  Viscoelastic  problem  reduces  to  an  associate 
elastic  problem  with  material  constants  as  functions  of  the  space 
variables . 

(c)  Temperature  uniform  but  time  dependent:  L.T.  applicable  with 

respect  to  the  reduced  variable  ^  Viscoelastic  problem 

reduces  to  an  associate  elastic  problem  as  in  case  (a). 

(d)  Temperature  both  space  and  time  dependent:  L  T.  inapplicable. 
Associate  elastic  problem  does  not  exist. 

A  certain  amount  of  simplification  is  achieved  in  cases  where  the 
dilatatlonal  response  of  the  material  is  elastic.  If  the  material,  is 
also  thermorheologically  simple,  then  the  dilatatlonal  response  remains 
unaffected  by  temperature  changes  (since  when  the  function  is  shifted 
it  merely  reproduces  itself). 

The  dilatatlonal  behavior  of  the  material  is  now  given  in  the 

fora 


s««  *  i  ©}  (54) 


where  K  is  now  independent  of  time  and  temperature . 
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An  Alternative  Integral  Fora  of  the  Constitutive  Equation* 

It  oust  be  emphasized  that  Eq.  (l)  and  (5*0  are  not  a  unique 
Integral  representation  of  the  viscoelastic  constitutive  lavs.  One. 
for  instance,  may  make  use  of  the  "tension  modulus"  E  and  "poissons 
ratio"  v>  ,  where  from  a  purely  formal  standpoint  E(t)  is  the  stress 
response  of  a  tensile  specimen  to  a  constant  axial  unit  strain,  aid 
V  ( t)  is  the  corresponding  lateral  strain  response  to  the  sane 
strain . 

Under  isothermal  conditions  the  stress-strain  relations  for 
direct  stresses  and  strains  become , 

-t 

"el*  ^  *6s«)  (55) 

a't  >o 

where  suffices  1,  2,  3  are  cyclically  interchangeable. 

For  a  material  with  elastic  dilatational  response  the  formal 
equivalence  of  (5*0  and  (55)  can  be  found  by  taking  L. T. 

Then  (l)  and  (51*)  become: 

5  -  =  J>  G  e- 

G  = 

K  K 

whereas  (55)  now  is: 


(56) 

(57) 
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(58) 


16 


or 


-  (I-+  Ob)Zu  _  (59) 


£ 

ii 


1+ 

\>  E 


(60) 


Comparison  of  (56),  (57)  and  (60)  and  a  few  algebraic  manipulations 
yield 


9  Kg 
Q>  *+-  |>  G 


(61) 


\7  _  _L  =SK-  I*  & 
*  GfT  +  f>G 

G  =  E  _ 

'-+  (>  v> 


(62) 


(63) 


K  r  }>  ^ 

3  (|  -  2|>^>) 


(64) 


Hence  in  terms  of  E  and  V  the  complete  stress-strain  relations 


are: 
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Since  all  properties  should  obey  the  sane  shift  Ism  (65)  and  (66) 
can  be  generalized  for  the  non- isothermal  case  l.e. 


d'C  - 


srnd  naturally 


(c,j>  (67) 


C  GCf-r)9-- 

Jo  T>  * 


cL*t 


(£*i> 


(68) 


In  the  subsequent  analysis  ve  shall  have  occasion  to  use  both 
constitutive  representations. 
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CHAPTER  I 

Thermal  Stresses  In  Viscoelastic  Solids  With  -Material  Properties 
Exhibiting  Thermorheologically  31mple  Temperature  Dependence. 

1.1  Introduction 

It  was  pointed  out  In  the  general  introduction  that  when  the 
temperature  field  is  transient  and  space  dependent ,  elimination  of  the 
time  dependence  from  the  constitutive  and  field  equations  simultaneously 
has  been  Impossible  by  means  of  L.T.  or  any  other  exact  transformation. 

Part  of  the  present  chapter  will  be  devoted  to  developing  approx¬ 
imate  techniques  by  which  time  dependence  may  be  eliminated,  thus  facil¬ 
itating  the  solution  of  the  relevant  equations.  In  fact  two  such 
techniques  are  given  which  are  not  dependent  on  the  particular  geometry 
of  the  body  at  hand. 

However,  other  approximations,  such  as  material  incompressibility 
or  effective  constant  Poisson's  ratio,  that  are  suited  to  a  particular 
geometry,  such  as  the  hollcw  cylinder,  are  also  used  wi+v  -Wantage . 

The  inapplicability  of  the  principle  of  superposition  -  except  in  a 
very  restricted  sense  -  robs  the  analysis  of  a  powerful  tool,  and 
limits  the  chances  of  development  of  an  exact  general  theory. 

1.2  General  Approximate  Techniques 

An  approximation  rendering  L.T.  applicable 

Consider  the  devlatoric  constitutive  relation  ( 35)  in  the 
presence  of  a  transient  non-uniform  temperature  field. 
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(1.2.1) 


when  the  suffix  is  emitted. 

Figure  (1.2.1)  shews  the  variation  of  ^  with  tine  for  aono- 
torically  Increasing  and  decreasing  temperature. 


Figure  (1.2.2) 
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The  present  technique  hinges  on  the  piecewise  linearisation  of 
over  saall  tine  intervals, 


Then  in  any  typical  Interval  K, 

?  =  A«  t  *  ,  *  J  t  <  +„  (12  2) 

whar.  Is  the  "elope"  of  £  in  the  Kth  Interval  (Pig.  1.2.2) 

We  new  as suae  that  by  the  present  procedure,  we  have  found  the 

expressions  for  *?-;•(£ )  for  each  interval  up  to  t*.  ,  and  call  the 
'  K 

expression  for  -€t  j  In  the  Interval  K  by  -€<■  (t)  where  t  S  ^  *  . 

By  introducing  (1.2.2)  and  carrying  out  the  integration  in  (1.2.1) 
over  all  intervals  for  (:  ^  «C  and  assuming  that  fVf)  continues 

as  a  straight  line  for  t*.  <  t  $  od  we  obtain: 

K. 

«  =  > 


(1.2.3) 


(1.2. *0 
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Then 

K. 

*;;(0  r  ^  I  (1-1  »  I  O) 


We  now  make  the  following  transformation 

x  - 


where  now  o  <  *,  <  <x J 

and  in  the  last  integral  T 

—  *-*- 1 

y  =  *-4*. 


where 


o  $  ^  C 

In  view  of  (1.2.7)  and  (1.2.8) 


€c.  (.*)  =  h  +  t«.)  = 


*+1 


K+l 


*■*-1 


*-4-l 


(oL*tV>  =  -£f;  (*-) 


(1.2.5) 


(1.2.6) 


(1.2.7) 


(1.2.8) 


(1.2.9) 


(1.2.9a) 


Ob  substitution  in  (1.2.4)  ve  obtain: 
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M+l 


Jo  2*1 


(1.2.10) 


How  taking  L.T.  of  (1.2.10)  and  letting 


i«.,  (h  -  L  l  <ii».,w}  <is-u> 


we  obtain 


J~,  ^  r  L,  QlL)fcu'- <i2“> 


But 


k+l 


i^=o)  =  ~e^U„)  ~  €:■  (1.2.13) 


in  view  of  the  continuity  of  strain  at  and  la  known 

from  the  previous  solution. 

Thus  (1.2.12)  be cooes 

-.jt° = L,  \ y  C- <■ <t j)  (i.8. 


1M 


The  expression  (1.2.4)  for  lv  with  the  change  of  variable  (1.2.7) 

IV 


be cooes : 
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since 


*  K-1 

K\a3  J* 

K  '9-* 

*< 

/  H 

ik-i 

(1.2.15) 

^X+l  ^  **•  ■+  ^  K+-  1  =■  ^  M. 

(1.2.16) 

*< 

r  •< 

(1.2.17) 

Taking  L.T.  of  ( 1.2.16)  with  respect  tn  x 

*  i< 


Ijf )  -  [  <^1  ^  G I V.1 -  4* 4  i  *■*'  «•* / <*> 

Jo  in-' 

**  Ml 

L  ia 


- 

J  V  <4 


^*1 — *  l*_  .  f  _  Aa  « 

.1  <  G(t  ) L4-'  »«.■  «>* 
A.*>  *"*'  ** 


I(f)=;  « .  G(. 

K  4--^.  A**i 


(1.2.18) 


Since  ©Cj  (t)  for  t  5  i  £  ^  ,<  is  known  frc«  a  previous  step 

of  the  solution,  T  is  known  function  of  p  found  fro«  Eq.  (1.2.17). 

K 


2k 


Thus  calling 


(1-2.19) 


and  taking  L-T.  of  (1.2.3)  and  in  view  of  (2. 3.1*0  and  (1.2. IB)  and 
(1.2.19)  we  obtain: 


L  K+\ 

L  Gr  (—  ) 

A—  '4...  ~  1 


(1.2.20) 


I 


Eq.  (1.2.20)  la  a  linear  relationship  between  the  transformed 
functions  St-(-  (|>)  and  [f>~)  .  If  (t^  is  known  then  S;j 

H-4-1 

is  also  known  and  thus  <  >  •  (V)  and  hence  <  £ :  can  be  found 

a*' 

from  (1.2.20)  by  an  inverse  L.T.  Therefore  the  strain  ■  (t)  is 
known  as  a  function  of  time  in  the  interval  t  £  V  ^  *<*  . 

This  solution,  however,  is  utilised  only  for  t  5  t  . 

For  t^(<t  the  curve  ^  (O  is  continued  with  a  new  slope  as  a 
straight  line  to  infinite,  and  the  above  procedure  is  repeated. 

Eq.  (1.2.20)  can  be  used  in  conjunction  with  the  solution  of  the 
boundary  value  problem  in  the  following  manner. 

In  the  interval  £  oq  (1.2.20)  is  written  in  the  form: 


(1.2.21) 
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where  i»  *  known  function.  A  similar  expression  nay  be  derived 

for  the  hydrostatic  constitutive  relation,  l.e. 


The  field  equations  for  the  current  interval  aay  also  be  put  in 
the  fora 


-4  44  ; 

•  \ 

\  Si 

-V  *  9 

1  J 

=  T  • 

O  M> 

/V  t 

-  (O 

* 

(1.2.23) 


Eq.  (1.2.21),  (1.2.22)  and  (1.2.23)  are  necessary  and  sufficient  for 
the  solution  of  the  non- isothermal  viscoelastic  boundary  value  problem 
in  the  transformed  plane. 

This  technique  has  been  used  with  encouraging  results  in  the  case 
of  the  slab,  see  Ref.  [.lo]  .  In  the  same  reference,  the  problem  has 
been  formulated  where  the  constitutive  equations  are  given  in  the  fora 
of  differential  operators. 


Reduction  of  the  Constitutive  Integral  Equations  to  a  get  of 
Simultaneous  Algebraic  Equations 


Consider  the  following  constitutive  equation  of  the  Volterra  type 


where  s  it  the  deviatoric  stress  tensor  and  e  is  the  deviatoric  strain 


tensor. 

For  simplicity  denote 

5  (1.2.25) 

then  (1.2-24)  becomes 

t 

(  G(*K,t,  *)*£  (1.2.26) 

J  O 

For  fixed,  €(0  can  be  approximated  by  a  piece-vise  linear 

function  of  t  (Figure  2.1.3)* 


Figure  (l  -t*3) 

LINEAR  PIECEWISE  VARIATION 
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Let  A  K  represent  the  "elope"  in  the  n^  interval; 
the  values  of  %  and  -€  at  the  end  of  the 


Then, 


Let 


-U 

C  G  (*  1 cW  =  G 


where 


4 


•V-l 


Then 


G 

M  V 


a 


where 


a  -v  =  A  „  A  ^ 


or  in  aatrix  fora 


i 


S 


[GJ  [4 


and 

interval. 


(1.2.27) 


(1.2.20) 


(1-2.29) 


(1.2.30) 


(1.2.31) 


(12.32) 
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where 


L<0 


Gv 

&',*>  6.,%l  G,tU 


L^Mu 


Note  that  |  6>J  i«  a  triangular  aatrix. 
Since 

a  =  -e  -  e 


we  have  the  relation 


=  CHJ 


where 


[Hj  . 


I 

I  I 
I  I  * 


I  I  1 


Fro*  (1.2.32) 


M  = 


and 


e  ^  LHJLGJ 


(1.2.33) 


(1.2.3M 


(1.2.35) 


(1.2.36) 


(1.2.37) 


(1.2.38) 
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•Iso 


H  •  [g.HhJ-'H 


(1.2.39) 


[?J  being  triangular,  J  is  very  easy  to  find.  Eq.  (1.2. 38)  is 

a  numerical  solution  of  the  Volterra  equation  (2. 3-2U) .  However,  more 
will  be  said  about  this,  in  the  closing  parts  of  this  chapter. 


Application  to  the  Hon- Isothermal  Viscoelastic  Boundary  Value 
Problem 

Equation  (1.2.32)  can  serve  as  a  basis  for  solving  the  non-isother- 
mal  boundary  value  problem  for  all  t  by  solving  consequentively  for 
small  finite  time  intervals. 

In  this  way  the  time  dependence  is  eliminated  and  the  problem 
reduces  to  the  solution  of  a  related  elastic  problem  with  "initial 
strains"  and  spatially  dependent  elastic  properties. 

In  the  general  problem,  parallel  to  (1.2.32)  there  will  be  a 
dilatational  stress-strain  relation  vhlch  again  can  be  put  in  the  form 


fO  s0\|[hH  il8-uo) 

also 


ve  call 


(gj  c«  r 1  r  lg  *j 


(1.2A2) 
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and 


(1.2.43) 


Where  j  J  and  £  Gj 


are  triangular  matrices. 

In  the  first  Interval  the  stress- strain  relations  becc 


(1.2.44) 


& 


I 


(1-2.45) 


Clearly  (1.2.44)  and  (1.2.45)  together  vith  the  equilibrium 
equations,  strain  displacement  relation,  and  boundary  conditions  consti¬ 
tute  an  elastic  problem  with  spatially  dependent  elastic  constants. 

In  the  second  interval  ve  get 


G*  ■€ 

l, Xl 


+  G 


(1.2.46) 


C*=  (1.2.47) 

We  new  have  an  elastic  problea  vlth  "initial"  Known  strains 
€,  ,  &  (  and  0,  .  Therefore  in  any  typical  interval  we  can  solve 
for  the  stresses  in  that  Interval  in  terms  of  the  temperature  in  that 
same  Interval,  and  the  stresses,  strain  &  temperature  in  all  previous 
intervals.  Itoe  solution  for  all  time  can  be  determined  in  this 
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systematic  way. 

It  Is  interesting  to  note  that  wring  boundaries  can  be  dealt  vlth 
by  aeans  of  the  above  wthod,  in  the  light  of  saall  deforwtion  theory. 
This  is  simply  done  by  dete raining  the  nev  boundary  froa  the  solution 
in  the  previous  interval  and  asauaing  the  boundary  to  rewin  fixed  in 
the  current  interval.  This  naturally  ia  only  approximately  true. 

The  above  approach  provides  a  systewtic  wthod  through  which  the 
general  viscoelastic  boundary  value  problem  can  be  solved  by  numerical 
wans. 

A  Variation  on  the  Previous  Technique 

An  alternative  approach  ia  to  consider  as  an  escalator  function 
of  tiw.  See  Pig.  (1.2.4) 


Figure  (1.2.4) 


Then 


-ett) 


-  HO-O  (-»=  0,1,2...^ 


(1.2.46) 


Substituting  (1.2.49)  In  (1.2.26)  ve  obtain 
__  t 

S(M  -  (1.2.50) 

*  o  _ 


Thus 


V  =  ^9'6(»,.{J,)=£fi  &  (1.2. 

o  „  _  f  ** 


51) 


where 


6t 


-  G  (n. 


(1-2.52) 


Eq.  (1.2.51)  Is  Identical  in  fora  with  Eq.  (1.2. 30),  however,  the 
coefficients  Gi)Kr  are  now  defined  differently.  Hie  advantage  of 
this  particular  formulation,  at  the  expense  of  smaller  time  intervals, 
is  that 


G 


l,KK 


G>  O') 

I 


(1-2.52) 


for  all  t'v  and  all  temperature  variations,  this  simply  being  the 
elastic  response  of  the  material  which  the  shift  hypothesis  renders 
Independent  of  temperature. 
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Haturally ,  we  bare  again  in  aatrix  fora: 

[<?}  (1.2.53) 

vhere  now  the  element*  of  the  coluan  matrices  i  5  ^  and  { -C  j  are 

the  average  values  of  the  deviatoric  streaa  and  strain*  in  the  corres¬ 
ponding  tiae  Intervals.  Also, 

{<2  }  =  [<**]  [h]  (1.2.54) 

Again,  in  the  first  ^tiae  interval  we  have  Eq's  (1.2.44)  and 
(1.2.45).  6tjl,  and  ()  ,  however,  are  now  constants  and  the 

viscoelastic  problem  la  reduced  to  an  elastic  problem. 

In  view  of  (1.2.52)  the  viscoelastic  problea  is  repeatedly  reduced 
to  an  elastic  problea  in  all  the  intervals  except  that  for  Intervals, 
subsequent  to  the  first,  known  "initial  strains"  will  also  be  present. 
The  power  of  this  technique  does  not  need  emphasising. 


1.3  Solution  By  Assuming  Mechanical  Incoapresslblllty 


This  assumption  has  been  utilized  in  the  past  2  to  obtain 
solutions  to  viscoelastic  boundary  value  problems .  It  is  based  on 
the  hypothesis  that,  where  the  dllatatlonal  response  is  elastic,  the 
bulk  modulus  Y\  is  of  exceedingly  high  order  of  magnitude  so  that 


in  the  expression 

£  =  —  +  3  *e  © 
3K 


(1.3.1) 


< 

! 


3* 


the 


<o 


ter*  is  negligible  end  hence 

e  =  3  O/a  0 


(1-3.2) 


This  simplification  can  lead  in  certain  caaes  to  a  closed  for* 
solution  -  see  sections  h  and  5a  -  however,  unless  ^  is  in  fact 
very  large  such  a  solution  can  only  be  approximate. 

On  the  other  hand,  such  a  solution  can  serve  a  very  good  zero' th 
approximation  to  an  iteration  process. 

For  the  sake  of  illustration  ve  give  an  account  of  this  process 
for  a  plane  strain  problem  ,  however,  the  method  is  quite 

general.  A  solution  is  found  for  ^  -  O  .  Fro*  this  solution, 

<o  is  calculated  fro* 

t  * 

=  £x-*-£-y  -4  ^  -  °o  (l«3*3) 


This  value  of  (£  is  substituted  in  (1.3.1)  and  I'C  is  taken  as 
the  actual  modulus  of  the  material.  Then  ©  is  found  in  terms  of 
a  new  temperature  where: 


©,  -  © 


j. 


(1.3.M 


where  £>0  is  the  v*lue  found  fro*  the  first  solution  and 


t,  = 


(1.3*5) 
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A  second  solution  is  now  found  in  terns  of 
repeated  to  required  accuracy. 


and  the  process  is 


1.4  Application  to  a  Viscoelastic  Cylinder  Enclosed  in  an  Elastic  case 
Formulation  of  the  Problen 

Here  we  consider  a  hollow  cylinder  under  plane  strain  subject 
to  a  radial  transient  tenperature  field.  The  cylinder  is  enclosed  by 
a  thin  elastic  shell  which  is  assuned  to  be  rigidly  bonded  to  the 
cylinder  (Fig.  1.4.1). 

Under  these  conditions 


I  '  T(rr,  t) 


(l.fc.l) 


M  -  4X  (<1  t)  (1-4.2) 

or  » 


(1-4.3) 


9u. 

'd'f 


(1-4.4) 


We  now  have  the  following  compatibility  relations  in  view  of  (1.4.2), 
(1.4.3)  and  (1.4.4). 


£ 


-v 


9 


£ 


=  1  ? 
v  'i-y 


(1-4.5) 


(1.4.6) 
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£•»•  -  _  £-»  -  €q  'b  €© 


*dv 


,  6=6  -«■£ 

■>  -*  © 


(1.4.7) 


where  and  -€e  are  deviatorlc  strains. 

The  equilibrium  relation  to  he  satisfied  is 


6y-6© 


-  O 


(1.4.6) 


(54J 

Making  use  of  Eq.  we  get 


t  *  -  ■* 

i  K 


(1.4.9) 


Iteration  Solution 

We  new  empl9y  the  hypothesis  that  the  material  is  elastically 
incompressible  in  dilatation  by  mechanical  forces,  £.2^  ,  i.e. :  • 


£  _  o  (1.4.10) 

and  use  this  assumption  as  the  zero' th  approximation  to  the  solution 
of  our  problem,  as  in  the  case  of  the  infinite  slab. 

Then  from  (1.4.6)  and  (1.4.9)  and  (1.4.10): 


6 


'i  w 


SK' 


(1.4.11) 
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Integrating 


9 


-Sir 


(1.4.12) 


Making  uae  of  (1.4.10): 


(1-4.13) 


or 


(1.4.14) 


where 


(1.4.15) 


Fro«  Eq.  (1.4.7)  we  obtain: 


34e 

2^ 


(1.4.16) 
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which  In  view  of  (1.4.9)  and  (1.4. 10)  becoaes; 

-  i  ^^o©-2.60V  (1.4.17) 

jv  -v  l  J 


or 


'th  d-"- 


18) 


We  now  uae  (1.4.7),  (1.4.8)  and  (30)  to  obtain 


-t 

- „f G($-?')  2( -'V* 


(1.4.19) 


which  in  view  of  ( 1.4.18)  becoaea 


3V 


!*(  G(f-f') 

-»  J0  2* 


w’  )o 


a-t 


ak't 


(1.4.20) 


-  —  1  G(f ')^  (2'Vjj_©)d^  (1.4.21) 
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Then 


t 

*.1  (&$-$') 

->*  K  a* 


(1.4.22) 


and 


"t  -r 


o 

I 


3-V 


d>  dt 


(1.4.23) 


Obviously  -  q  since  or-^  ,  is  a  free  surface. 

If  at  -y  there  is  a  rigid  boundary  £  (vt  ,t")  -O  and  the 

complete  solution  is 


£ 


Kp,  0  c\j> 


(1.4.24) 


In  the  case  of  the  boundary  being  an  elastic  shell  6-^  and 
for  both  cylinder  and  shell  must  be  equal  at  the  conaon  boundary. 

Let  (S?s  (t)  be  the  temperature  of  the  shell  assumed  constant 
through  the  thickness.  All  parameters  with  suffix  “5  refer  to  the 
elastic  shell. 

Then 


^©K,tr)  =  0+O«V«  ©s  - 


(1.4.25) 
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Also  let 


*  X 

— 1  £»f f ^  =  G(t,fc) 


(1.4.26) 


Then  (1.4.23)  becomes: 


F  (*,  *■  )  d 


*r  -4 


-t 


2* 


(14.27) 


Eliminating  &nt  from  (1-4.25)  and  (1.4.27)  we  get 

F  , 

E<  ^  U.  4-  (  t-Cr,  OdT  r  fe. _K,0 

0-0^  ® 


(1.4.26) 


Eq.  (1.4.28)  is  a  Volterra  Integral  equation  of  the  second  kind. 

The  solution  of  this  equation  -  see  Appendix  II  -  for 
and  substitution  for  0  in  (1.4.23)  yields  the  complete  solution 

of  the  problem  in  the  zero'th  approximation. 


The  circumferential  stress  de  is  obtainable  from  ( 1.4.8)  i.e. 


(1-4.29) 


To  determine  (£  we  make  use  of  the  identity 

4=3  (<U-0 


( 1-4-30) 


where 


(l-1*. 31) 


Now 


e  ) 


-  * 


(1-4-3 2) 


In  view  of  (1.4.30),  (1.4.31),  (1.4.32)  and  (1.4.23),  &0  can  be 
determined. 

Substituting  for  d>0  in  (1.4.9)  and  repeating  the  above  proce¬ 
dure,  we  can  find  a  first  approximation  tothe  solution  in  tens  of 
some  fictition  temperature  0,  where 


A  systematic  formulation  of  the  Iteration  process,  together  with 
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sufficient  conditions  for  its  convergence  is  given  in  Appendix  II. 

A  computer  progrsn  based  on  this  theory  has  been  developed  9 

{•uadi 

by  aeans  of  which,  the  stress  and  strain  distributions  are^in  a  hollow 
cylinder  contained  in  a  rigid  shell.  The  results  and  the  corresponding 
tenperature  histories  are  shown  in  Pigs.  (1.4.2),  (1.4.3),  (1.4.4). 
(1.4.5),  and  (1.4.6). 


TIME(HR) 
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HOOP  3TRS3S  VERSUS  TIME  FOR  STEP,  STEP  WITH  A  JOG, 
AND  CYCLIC  INPUTS 
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1.5  A  Better  Zero’th  Approxlwatlon  to  the  Solution  "ab  lnltlo" 

By  a  slightly  different  approach  a  closed  fora  solution  is  obtained 
vhich  involves  a  less  drastic  approxioatlon  than  the  one  eoployed  in 
the  previous  section. 

Consider  the  following  elastic  relations  vhich  can  be  lmediately 
generalized  to  the  corresponding  viscoelastic  relations. 

r  £  \  ~  °  (  **  +  6y)}  +  **  ©  U-5.1) 

Since  G.  -Z.O  we  get 

-  E«/0  ©  (1.5.2) 

and  hence : 

d  =•  (l-<- ^)  (  d-,+6e>)  -  Efl/„©  (1.5.3) 


However 


(l  V  ) 


(1.5A) 


and  the  corresponding  viscoelastic  relation  for  constant 


K 


is 
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Since  HC  is  conetant  wad  E(t')  varies  between  Eo  and  tero- 
for  uncrossed  linked,  and  E Q  and  E- ^  where  E^  2;  ^  E.^  for  a 
real  propellant, 

,  +  v  ^  i+wo  <:  i  (1.5.6) 

2 

For  polnethyl  methacrylate  v>  =  .35,  which  is  a  low  value  for 

© 

viscoelastic  materials.  However  even  in  this  case 

l-SS-  <  I  4-  O  ((-■)  ^  l  -S  (1.5.7) 

i.e.  varies  within  narrow  limits  and  thus  one  is  Justified  to 

make  the  approximation 

I4-n?(0£  ELLLL5'  -  l-4-2.r  u?  (1.5  8) 

z 

where  ^  «  .425  in  this  case. 

Then  by  considering  a  constant  average  v?  for  a  viscoelastic 
material  we  will  have  an  expression  analogous  to  (1.5*3): 

4  ,  (1.5.9, 

'*>  o'*- 


In  view  of  (1.4.8)  we  get 
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6  + 


or 


to  ~  { I  ■ 


also 


r  >  '-  h  -  ,0<°® 


Combining  (I.5.H)  and  (1.5.12) 


i  -  ("\)  =  —  -  -G‘<S  )-!i:(f(f-f'!S®dii 

-*a*  sl  3k  ^j*v  v> 


Integrating  (1.5. 13) 


Vfe-’  £R'’‘M£r(E(*-*,)!?<“ 


j. 


3* 


|>  ©  dp  -+  C  (0 


(1-5.10) 


(1-5-11) 


(1-5-12) 


(1-5-13) 

H© 


cr 


& 


(1-5-14) 
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and 


In  view  of  (I.5.I6)  and  (l.U.15)  we  get: 

t  r  ]±Z  £  -  —  X2  +  *OP  +  ^ 

0  **  =}K  V 

where  _Q  (^")  =  =  O 

We  now  differentiate  (1.5-17)  with  respect  to  •f’  and  denote  such 
differentiation  by  a  dot,  then  after  differentiation  with  respect  to 
tiae,  Multiplication  by  Gi  )  ,  and  integration  between  the 

Halts  of  o  and  ~t  ve  get: 


(1-5-17) 

(1-5-18) 
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by  (1.4.19). 

Let  us  consider  now  the  functions: 

?  =  §  fa  O  ,  ^  =  f  G*.  o 

For  fixed  -v  ,  ^  Is  a  aonotonic  function  of 

It  can  be  inverted  In  the  fora 

'  <■>  iKf) 

Then  if  ^  Is  a  function  of  or  and  "t  , 

f  -  {b,  =  fK£ ) 


d* 

0®^* 


(1.519) 


(1-5.20) 

and  hence 


(1-5.21) 


(15.22) 


Also  let 
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-  =  Zko  -  7k?) 

^  —  - 

^5  =*  ct>K  O  =  cJ?M  ) 

9'v 

^  s  TTtv,e)  =  1TK?) 

Then  substituting  (1.5*23)  in  (1.5*19)  ve  obtain 


4  4 


(1.5.24) 


Referring  all  quantities  to  the  ^  plane  ( 1.5*24)  be coaes 


(1*5*25) 


We  new  take  Laplace  transform  of  ( 1. 5.25)  with  respect  to 


and  obtain: 


(Vo 

=5^ 


i»  & 


—  j\ 

-+-  3  G  1 1 


_2 

v3 


fi  K, 


Now  letting 


(1.5-26)  becomes : 


(1-5-26) 


(1-5-27) 


(p-3c/mf>#  TT  +  3-  ^  C 


(1-5-28) 


Now  taking  the  inverse  Laplace  transform  ve  find 


5  6 


Finally  reverting  to  the  t  plane  (1-5-29)  take*  tha  Tom. 


+  1  f  (1‘5-30) 

'V3  '0  3* 

Substituting  for  £  t  d> ,  TT  from  (1.5-23)  and  after  putting 


(  -  3  J- ot't  ^.O  d.5.31) 

•  A 


we  obtain 


(1-5-32) 


Hence : 


Sde  di 


(1.5.33) 
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since  £  nr  =  being  a  free  surface.  Therefore, 

<UvO  +  (i.s.» 

%  or,  o 

Also  in  view  of  (1.5.17)  and  (1.5.18), 


fen  =  a  -+  ?±Le>)  (1.5.35) 

®*  3K  ^ 


=^eKvO  5  6^  =  £  (-v  ,  f) 


(1.5.36) 


At  the  interface  of  the  viscoelastic  cylinder  and  the  elastic 
shell  we  have  the  relation 


0^VS  (S>s  -6 


^  k 


(1.5.37) 


where  suffix  5  refers  to  the  shell  and  is  the  nean  temperature 

over  the  shell  thickness. 

Eliainating  from  (1.5-35)  and  (1.5.37)  we  obtain: 


i±*  4  .+  cj1^  -  -t  lz?a }  (1-5>38) 

3*  ***  */  Es 
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In  view  of  (1.5.3**)  and  (1.5*38) 


-2 

Y 

1 


$  I +)? 


1-3 


1  ^  v  c(H-^  r  ?£' 

u.  J  1  ) 


(1.5.39) 


where 

"a 

S^Ct.t)  r  f~  Rf?KO  (1.5A0 

Eq.  (1.5.39)  i«  again  a  Volterra  Integral  equation  of  the  second 

kind  in  C?  (0  . 

1 

Solution  of  this  equation  and  substitution  in  (1.5. 33)  solves  the 
problem  completely. 

The  hoop  stress  can  be  found  immediately  from, 


£ 

© 


£  +  ^ 
^  l>r 


(1.5.41) 
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Since  this  solution  is  dependent  on  the  aforeaentioned  approxima¬ 
tion  one  can  test  its  accuracy  by  calculating  id  from  the  expres¬ 
sion^ 


and  then  use  the  iterative  method  developed  in  this  chapter  to  calculate 
a  nore  exact  value  of  .  However,  comparison  of  6  as 

obtained  from  (1.5.42)  as  opposed  to  (o  obtained  from  the  approxi¬ 
mate  expression  (1.5*9)  will  determine  whether  an  iteration  should  be 


necessary. 
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CHAPTER  II 

Dynamic  Stresses  la  Thermorheologically  Single  Viscoelastic 
Bodies 

2.1  Introduction 

Dynamic  stresses  in  linear  viscoelastic  elastic  solids  under  non- 
isotbermal  conditions  are  still  an  unexplored  field.  Hopes  of  obtain¬ 
ing  closed  form  solutions  even  for  the  simplest  configurations  are 
rather  small  as  will  be  appreciated  from  the  contents  of  this  report. 

One  will  recall  that  the  quasistatic  problem  of  the  slab  and  the 
sphere  were  given  a  closed  form  solution,  formally  at  least  because 
the  relevant  equilibrium  equation,  in  both  cases,  could  be  integrated 
directly.  When,  however  inertia  forces  are  taken  into  account  this 
is  no  longer  possible. 

In  the  present  chapter  progress  has  been  made  by  limiting  our  atten¬ 
tion  to  incompressible  viscoelastic  materials,  in  the  sense  that  volu¬ 
metric  changes  either  due  to  mechanical  forces  or  temperature  fields 
are  zero. 

Consequently  the  dynamic  stresses  examined  here,  arise  because  of 
the  time  wise  variation  of  the  mechanical  forces  (stresses)  applied  at 
the  boundary. 

We  limit  ourselves  to  the  configurations  of  the  sphere  and  the 
infinite  hollow  cylinder,  both  with  polar  symmetry,  so  that  dependent 
variables  are  functions  of  the  radius  only  and  time.  On  the  other 
hand,  within  this  restriction,  temperature  fields  are  both  non- homo¬ 
geneous  and  transient  in  nature. 
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The  solutions  of  both  problems  reduce  to  Volterra  integral  equations 
of  the  second  kind,  which  can  be  systematically  solved  numerically  with¬ 
out  undue  difficulty. 

In  the  first  part  of  the  report  we  treat  the  viscoelastic  sphere 
with  polar  symaetry.  In  the  second  part  we  treat  the  hollow  viscoelastic 
cylinder  infinite  in  length  and  also  with  polar  s yaw  try.  In  both  cases 
the  dynamic  stresses  are  exaained  in  the  presence  of  ties- varying  mechan¬ 
ical  forces 

2.2  Dynamic  Stresses  in  a  Viscoelastic  Sphere 

Sphere  With  Infinite  External  Radius 

We  first  Investigate  the  case  of  a  sphere  with  an  infinite 
external  radius,  so  that,  essentially,  we  have  a  viscoelastic  continuum 
with  a  spherical  cavity,  in  a  transient  inhomogeneous  temperature  field 
with  polar  symmetry.  The  cavity  is  subjected  to  mechanical  pressure 
which  may  vary  with  time.  It  is  assumed  that  prior  to  time  t-0  the 
pressure  was  constant  and  that  any  variation  began  at  time  t  ■  0+. 

The  temperature  field,  however,  may  have  existed  prior  to  the 
application  of  the  pressure,  the  time  of  application  of  the  temperature 
field,  relative  to  that  of  the  applied  pressure  variation,  being  arbi¬ 
trary. 


Figure  (2.2.1) 
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Let  ,  be  the  pressure  In  the  cavity,  "a"  the  radius  of  the 

cavity,  -f  the  radial  distsnce  froa  the  center  of  the  cavity,  and 

U.  the  displaceaent  along  w 

Then  in  the  usual  notation,  and  in  view  of  the  assumed  polar 
syaaetry  ve  have, 


£  .  H 

'b* 


6  =  £  -  H 

'  e  4>  < 


(2.2.1) 


6=6  ^  .  ZU.  K  = 

<  K  '  +  -  ) 

3-v  -v 


(2.2.2) 


(2.2.3) 


6  1.  -  o 

K 


(2.2.4) 


The  equilibrium  relation  is 


**7  -  KV<V>  =  p~\ 

'Jv  "V  ‘  dl 


1>\ 


(2.2.5) 


Froa  (2.1.4)  and  (2.2.2)  ve  get 


2a  2a  _ 
-r 


(2.2.6) 


V-  -  C(0  ^ 


- 1 


Hence 


(2.2.7) 


where  Q{t)  la  a  constant  of  integration. 
Also  in  view  of  (8.2.1) 
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£ 

t 


0 


-  -  l  C(t)  ■* 3  ) 
=  C  (O  *  / 


t  .  t  .  CCri 

'T  0 


Hence  froa  (8.2.3) 


6-6  = 
->  9 


dM: 

A 


(2.2.8) 


(2.2.9) 


(2.2.10) 


Substituting  in  (2.2.5)  we  get 
t 


at  + 1  *le 


(2.2.11) 


Integrating  (2.2.11)  and  in  view  of  the  fact  that  (oji*)  ■  -  j5^ 


we  get 


6  •+■  J>  it)  -  Co  j"  n'  ^  ^  G  ^  ^ 


«•  ^ 

Interchanging  the  order  of  integration  in  (2.2.12) 

6  +  j>(ri  =  (,  |  |f  -f  Wd 

+  *lc 

a  ->  j-£*- 


(2.2.12) 


(2.2.13) 
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Mow 


(2.2.14) 


provided  that  the  lap roper  integrals  exist,  and  they  do,  since  G  is 
-4 

bounded  and  V  is  saonotonically  decreasing  function?  Let, 


and 


(  ^/4 

Substitute  (2.2.14),  (2.2.15)  and  (2.2.16)  in  (2.2.13)  we  get 

j>(?)  r  <b|  ^  K'ct.v)  -  kJc+.'v.-t')  y  ^1  -fc 

4-1(1--)^  (2.2.17) 

Cl  *  d-i*- 


(2.2.15) 


(2.2.16) 


Mote  that 

f\J  0,  *  ,  —5*  O 


(2.2.18) 


as  *  <*0  . 

Also  since  as  ~t  -*•  oO  ,  (2.2.17)  becoaes 


*  At  rate  greater  than  -i  . 
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Q> 


i 


O 


cJ* 


_r  is 

o-  di" 


(2.2.19) 


Eq.  (2.2.19)  completely  determines  C  (0  ,  hence,  formally  at 
least,  the  problem  is  solved. 

Also  In  view  of  (2.2.19),  eq.  (2.2.1?)  becomes 


or 


6  r  - 


a 


%  - 


(2.2.20) 


C. 

(2-2-sl> 

o 

Hence  the  solution  to  the  problem  hinges  on  the  solution  of  (2.2.19), 


Let 


Then 


(2.2.22) 


(2.2.23) 
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aince  frc*  (2.2.7),  by  virtue  of  u.-*o  57=0, 

*  0% 


at  -t  =  o-t- 


c 


=  o 


-  o 


Then 


v 

=  (W,*) 


J, 


Je 


Hovever  regarding  (2.2.23)  as  a  double  integral  ve  get 


Let 


* 

<(*,*)  7 


'-Mf 


Figure  (2.2.2) 


r 


(2.2.24) 


(2.2.25) 


(2.2.26) 


<ct(  t)  d't 


(2.2.27) 
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Then 


DC 


(  Kd  %)  —  cM  = 
)  dt 


[  Kii,  cU'*<* 

Jm  )o 

■i  * 

K(t  (2.2.26) 


where 


Substituting  (2.2.22)  end  (2.2.26)  in  (2.2.19)  we  get 


|>(0  -  ^ ^  K  ct'0<p)<n  + 

1  o 

■i 

K\o  -  ^ 


(2.2.29) 


(2.2.29*) 


Equation  (2.2.29)  is  a  classical  Volterra  Integral  Equation  of  the 
second  kind  which  can  be  solved  numerically  by  various  Means. 


2.3  Temperature  Uniform  But  Time- Dependent 

In  this  case  eq.  (2.2.11)  can  be  directly  integrated  with  respect 
to  ~c  ,  to  give, 


+  -f  (<-*)  % 
a  n  D\ 


(2.^1) 


Since  &  -vo  as  ~t -v  vs  get 

•4 

i 


Mtt  r  1  (  Q  *  I  £C 

r  a?  )  «■  ^ 


(2.3.2) 
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or 


-t 


'4.?  ^  2-  ')d—  ^  =  *  b(e) 

Jl1-  aV)  <»*  r  1 

*  0 

(2*3*3) 

Let  ■* 

(2.3-M 

2  G»0 

^7  = 

(2*3*5) 

1 

where  Gq  •  0(o). 

Then  (2.3.3)  becomes 

•t 

(2.3.6) 

Also  fro*  (2.3*1)  end  (2.2.22) 

1 

-t 

(2.3.7) 

By  virtue  of  (2.3*2)  end  (2.3*7) 

(2.3.6) 

where  g(t)  is  given  fro*  (2.3.6). 
Fro*  (2.2.5)  end  (2.2.7) 
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2Ss 

-  f^)  V,  -  ^  i te)  C 1  *0‘J  <8'3-,) 


2.4  Sphere  With  Finite  External  Radius 

a.  Condition  of  Fixity  at  External  Radius,  -r  »  b 

If  u  la  to  vanish  at  r  «  b,  c(t)  £  0  frost  (2.2.7). 

Hence  u  =  0  and  all  deviatoric  strains  and  stresses  vanish  identically, 
The  equilibrium  relation  (2.2.5)  becomes 


(2.4.1) 


or 


(2.4.2) 


in  view  of  the  boundary  condition  it  r  >  a. 
Also  since 


(2.4.3) 


(2.4.4) 
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3tresses  are  transmitted  lnstaneously  and  there  are  no 
effect* . 

2.5  Free  Surface  at  -v  *  b 

Again  we  quote  (2.2.3),  i.e. 

i  ^ » =  4  £  S  « -*  '>  *''•*  +  b (,-^S 

'J  O  - 


at  Vob,  6*0  hence 


=  c{4* 


dA: 


Let 


^  G(?-f')'v7V/  =  L(4,V> 


Then  (2.5.2)  becomes 

i 


p(0 


(2.5.4)  is  an  equation  similar  to  (2.2.19). 

On  substitution  of  (2.5.4)  in  (2.5. 1)  we  find  for 

ih 

^  &  f  «!5  [  f  V*G(f-f')d^d* 

f  djc 

~  *  d-t* 


inertia 


(2.5.1) 

(2.5.2) 

(2-5.3) 

(2.5.4) 


(2.5.5) 


71 


2.6  Isothermal  Sphere  With  Infinite  Radius 
From  equation  (2.3.3) 


t 

4^  _L  (  d*  =  2- |»w 

at1  c-y  [  <**■  r 1 


Consider  G(t)  to  be  given  by  a  Maxwell  model,  i.e. 

£ 

GOO  :  fi9  ^ 


(2.6.1) 


(2.6.2) 


Taking  Laplace  transform  of  (2.6.1)  and  in  view  of  (2.2.24) 


s’c  4  -  Si  C  -  —  t> (sV 

*xr  r r 


From  (2.6.2) 


GiQs') 


substituting  in  (2.6.3)  ve  obtain 


(2.6.3) 


(2.6.4) 


S*C 

X.  ^  1 

a. 

'  r 

j>(^ 

(2-6.5) 

ar  Sj_5 

Let  S XC 

-  ^  (s)  and 

2.6?© 

aT 

-  Ul2 

(2.6.6) 

then  (2.6.5)  becomes 
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1 


CO  + 


-  £  bcs) 


(2.6.7) 


or 


v  “  .1C 


Sl  +  Y  +  U)1, 


(2.6.8) 


or 


<J  (O  -  -  f’Cs')  -  f,1**  jo  (0  • 

f  *■£♦“** 

Substituting  (2.6.9)  in  (2-3.8)  we  find 


(2.6.9) 


= 


(2.6.10) 


£ 


to 


U3 


X  -a* 


/  x  T“~ 

^XU>  -  — r 

4> 


e  />•»«. (2.6.u) 
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Then  (2.6.10)  becoaKS 

£  _  +  (2.6.12) 


vhlcb  is  the  complete  solution  for  all  load  cycles 


f(t) 


Zjo* 

r 


Figure  (2.6.2) 
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2.7  Dynamic  Stresses  in  an  Infinite  Hollow  Cylinder 

We  finally  consider  the  dynamic  stresses  arising  In  sn  Infinite 
cylinder  In  radial  transient  teaperature  fldld.  As  stated  In  the  Intro¬ 
duction  these  are  considered  the  result  of  varying  interval  pressure. 

Under  the  above  conditions,  and  in  vlev  of  the  assuaed  material 
lncoapressibility,  ve  have  in  the  usual  notation: 

e  =  “  t  z.o  (2.7.1) 

*»  h'i  ’  ®  -f  ’  4 

t  r  (  r  -  rO  (2.7.2) 

0  g-v  ^ 


Froa  (2.7.2) 


where  C(t)  is  a  constant  of  integration. 


Hence : 


b  _  3m.  =  _  1  C(H 

Tr 

^  ±  ^ 

(2.7.5) 

and 

fc-<=  =  ^  c«-) 

•v  ©  /y  *• 

The  equilibrium  equation  to  be  satisfied  is 

(2.7.6) 

'Dt-,  £-»  -6e  _  u- 

^  -V  1  9i' 

In  view  of  the  constitutive  relation 

-t 

(2.7-7) 

Sc:  r  (  Gd-'t)  ol* 

4  J. 

(2.7.8) 

and  the  the reo rheologic ally  simple  nature  of  the  cylinder 
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Substituting  for  4  froa  (2.7*6)  In  (2.7.9)  v«  obtain 


(2.7.10) 


As  a  consequence  of  (2.7*10)  end  (2.7-3) »  (2.7.7)  becoaes: 


^  -  GCf-V)  +  1 


*c 


w  ?  Jo 


7>\ 


f  9c 


(2.7.11) 


Free  external  surface 

Let  the  Interval  radius  be  "a",  and  the  external  radius  be  "b". 
Consider  the  following  boundary  conditions  <0  -  —  fKt')  at  <c  »  a  , 

V  1 

(,  *0  at  'f  -  b  (2.7.12) 

-V 

In  view  of  (2.7.12)  integrating  (2. 7- 11)  we  find 

Since  at  'i  ■  b,  (J  ■  0,  we  hare  froa  (2.7.13). 


(>Ct)  r  2^  d*  (2-7-HO 

~&C.  ac 

This  equation  completely  deteraines  — x  and  ^  ,  being  analogous 

to  (2.2.19) t  where 
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Mft.O  -  f  — , 

nr  ^ 


(2.7.15) 


Also  froa  (2.7.13): 


dri  dk 


■U  U-  SIS-!'1 

—  r  _i  d't 

J.  **  \ 

■f  sl‘  ^ 


(2.7.16) 


Incorporating  (2.7.14)  and  (2.7.15)  in  (2.7.16)  we  obtain 

4 

•J 

£,  *  -  2  *,'<')  ^  +  f^2  -?*»)  t  (2.7-17) 

where 


QC-.t.t)  =  (  G'H  ')cU 


U  *1 


(2.7.10) 


Cylinder  contained  in  thin  elastic  shell 

Let  the  suffix  s  denote  quantities  pertaining  to  the  shell. 

At  -f  *  b,  as  a  consequence  of  the  continuity  of  u  and  (£  at 
the  cylinder- shell  interface  we  hare: 


C(0 

b* 


(2.7.19) 
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Incorporating  this  condition  in  (2.7.16)  at  ■{  »  b  vs  find: 

i 

E‘  *  ©  -*-  *>«■)  r  £±.  (£  )  dt 

7  04  «  a  '  U  '  tv  \ 

l-v>s  '-^S  0  0  )o 


p  a  /  b  s  VC 


(2.7.20) 


Eq.  (2.7.20)  can  be  reduced  to  the  canonical  for*  of  the  Volterra 

Integral  equation  of  the  second  kind  by  the  following  procedure: 
dc 

'  jjice  c,  -rr  are  zero  at  t  *  o+,  we  have 


«  * 


C(0  = 


d*  It 


(2.7.21) 


C(0  r 


d  K. 

dt1 


(2.7.22) 


1  * 

(  H(fc.-t)  ^  d*  -  mV)  - 

\  a*  J. 

* 

where  *«,*')  d*' 

A 

Eq.  (2.7.20)  in  riew  of  (2.7*21)  and  (2.7*23)  becoaes: 

i 

^  ^S)  +)>fri  =  |{ 2 M «.t)*  ft-t)  &  >  -k  (;)'| 


(2.7.23) 


(2.7.2k) 
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where 

<j  (0  =  £  (2-7-25) 

Calling 

r  (2.7.26) 

,~v* 

v«  hare  fro*  (2.7*24): 

4 

p  X«J(^  )  ^j(t-)  ^  *  (*,  l)  <J  ft)0»^  ^  (>(04.^  4*®  (2.7.27) 

which  1*  the  canonical  font  of  a  Volte rra  Integral  equation  of  the 
second  kind. 

It  is  worth  noting  that  as  b  -»-*o  log  £  -*  «»  and  g(t)-*>>  o  , 
hence  u-*  0. 

As  a  consequence. 


^  *  60  =_(>(!-)  (2.7.26) 

l.e.  stress  wares  propagate  instantaneously,  without  any  subsequent 
vibrations  being  set  up.  A  rigid  external  shell  gives  rise  to  the  saw 


result 
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CHAPTTO  III 

Horizontal  Slump  of  a  Viscoelas tic  Hollow  Cylinder  Contained  In 
a  Thin  Elastic  Shell 


Synopsis 

In  this  chapter  the  problem  of  the  horizontal  slump  of  a  visco¬ 
elastic  cylinder  contained  in  a  thin  flexible  shell  resting  on  a  rigid 
plane,  is  solved,  the  solution  is  within  the  scope  of  small  deformations 
and  classical  shell  theory. 

The  viscoelastic  stress-strain  law  is  of  the  integral  type,  and  it 
is  used  in  its  most  general  form. 

The  solution  is  purely  formal  in  as  far  as  no  numerical  results 
are  given 

3.1  Introduction 

Ibe  problem  of  the  horizontal  slump  of  a  viscoelastic  hollow 
dylinder  contained  in  a  rigid  shell  has  been  solved  jjrJ  .  The  solution 
indicated  that  prohibitively  large  displacements  develop  for  large 
times.  It  is  reasonable  to  expect  that  a  flexible  shell  will  exert  a 
less  restraining  influence  on  the  cylinder  and  will  allow  therefore 
even  larger  displacements  to  develop. 

An  exact  formal  solution  tothe  problem  is  obtained,  within  the 
scope  of  small  deformations  and  classical  shell  theory,  in  the  form  of 
sine  and  cosine  series  in  the  angular  coordinate. 
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3.2  Analysis 

The  treatment  of  Isothermal  viscoelastic  boundary  value  problem 
is  facilitated  at  least  formally  by  removal  of  time  dependence  by  the 
application  of  Laplace  transform,  and  consequent  reduction  of  the 
viscoelastic  problem  to  an  associated  elastic  one.  As  a  first  step, 
therefore,  ve  obtain  the  solution  to  the  problem  of  a  hollow  elastic 
cylinder  contained  in  an  elastic  shell  resting  on  a  horizontal  rigid 
plane,  both  being  deformed  by  virtue  of  the  forces  exerted  by  their 
own  velgbt. 

We  further  assume  that  the  cylinder  is  "long"  so  that  plane  strain 
in  the  axial  direction  can  be  assumed. 

Under  the  above  assumptions  the  geoswtry  of  the  cylinder-shell  con¬ 
figuration  is  defined  by  Pig.  3.U.I. 

The  cosqplete  solution  is  obtained  by  solving  for  the  cylinder  and 
shell  separately,  and  satisfying  continuity  of  radial  and  shear  stresses 
as  well  as  radial  and  circumferential  displacements  at  the  cylinder- 
shell  interface. 


3-3)  Cylinder  Analysis 

Satisfaction  of  equilibrium  in  the  radial  and  circumferential 
directional  yields  the  following  relations  in  the  usual  notation  of 
polar  coordinates; 


^  I  'ir.e  +  f  ck  CjotO  -  0 


■V  ■&© 


(3.3.1) 


=  O 


(3-3.2) 
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K 

«C,-6 

(3.3.3*) 

s 

=  6 

(3.3.3b) 

e 

e 

X 

*  -tre 

(3.3.3c) 

vhere 


t 

0> 

+  6 

9  9 

} 

(3.3.3d) 

(3.3.1)  *nd  (3. 

3.2)  becoac 

h 

l  Zt 

"f-  —  - 

^  ^  ~s© 

bd 

4-  — 

+  PCI  CoS0  -  O 

(3.3.4) 

2>v 

-*  £© 

.-r 

2>^ 

r " 

'ase 

is 

a'k 

4-  or 

+  2t  4 

36 

b© 

fV*'*6  =  0 

(3.3.5). 

We  eliminate  the  dependence  on  ©  by  taking  finite  Fouries 
t rang fora.  Let 


tr 


'tC'y.e^tuM©  <30 

l-w 

ti 

(3.3.6a) 

£  W,  0)  Cot  1*0 

L„  " 

n 

(3.3.6b) 

r  -r 

r  ?  (1,0')  cx>s^©j© 

'-T) 

(3.3.6c) 
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TT 

-  n 

In  viev  of  the  symmetry  of  the  problem  S.,  ,  S&  ,  (o 
functions  of  ©  but  is  an  odd  function.  Thus  ve  have 
ditions : 


-  ^  =SK”')  =l 

dg  T>Q  &e 


-  Ko)  =  -  C 7,0  =° 
a© 


^Ko)  =  *  C".  *0  -° 


By  virtue  of  (3-3.6)  and  (3-3-7)  Kq-  (3-3-1*)  becomes: 


-i  ^ 

&y 


%  -  V 


—  - 


d<S 


r 

v 


o  ?  w  ^  I 


-"M 


and  Eq.  (3-3-5)  becomes 


2  t 


l*  © 


m  at  l 

n»  1 


(3-3-6d) 

are  even 

the  con-., 

(3-3-7a) 

(3-3-7b) 

(3-3-7c) 

(3-3-8) 

(3-3-9) 


We  nor  assume  that  the  material  is  incompressible  so  that  £  is 


04 


indeterminate  from  the  stress-strain  relatione.  Bence  eliminating 
<£>  from  (3.3.8)  and  ( 3*3*9) *  and  putting 

s  =  £  -£  *  s  (3.3.10) 

i0  •>  e 

we  obtain: 

K'v  41  +. ac  45  ^4^5  -o  (3.3.11) 

d^  ch” 

The  strain  displacement  relations  in  polar  coordinate  are: 

6  t  -  i ^  ^  1  3JT  y  _»  t  (3.3.12) 

-»  "  3-v  >  e  ^  ^  '  J-te'i  l v2©  w  3 

^  —  o  t  +  (  +  t  — o  t  ■+£  — O  (3* 3.13) 

*  >  -*  e  1  ^  e 

In  riew  of  the  material  incompressibility  the  stress  strain  relation 

is 

Sil  -  Sje;i  (3-3-1**) 

where 

r>  ,  .  f.  .  _  £•  WW 

‘  Y  (K  suaed)  (3.3.15) 


or  putting 
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**■  &ir 

=--2fees 


The  strain  compatibility  relation  is 


*-,1  (-re .'1=2^5) 

aV*  *-  »v*> 


or 


_  z 


&6  v 

iTV 


.*  2  &© 


=>  ew  -4-  2-r 


$\>1>  .aJ  | 

t  av>©  a©  * 


l! 

a© 


D~< 


New  let 


I 

t r 


u  -  — 


cr  = 


r 


-ia  (pr,  e')  to  s  «*  0  d  © 


-r, 

n 

L‘ 


O‘^*,0')  S<WM©d0 


Again  in  view  of  tha  existing  synwetry 


vfro)  =  ofr*)  —  O 


(3.3.16) 

(3.3.17) 

(3.3.18) 

(3.3.19) 

(3.3.20) 

(3.3.21) 

(3.3.22a) 


(3.3.22to) 
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Eqs.  (3-3-12)  through  (3-3-19)  la  the  transformed  plane  becos 

l 

*  'b'r 

5  =5  iii  ,  y  =  - 
>  e  -v  ^7  '  o  i  L  -v  d-*  ^  J 

(3-3-12a) 

L  -o 

2 

(3- 3- 13a) 

S  =  & 

(3.3.14a) 

^KK 

(3-3-15*) 

*  =  6  i 

(3-3-16*) 

S  *  lf,  6  =-lfo6 

(3-3-17*) 

cK  d-,**  cl'y  l  dv  d  J 

(3-3-19*) 

Substituting  for  S  froai  (3-3-16)  and  (3. 3. 17)  we  get 

4  4k  +  St\  =C  (3-3.23) 

dii  d^X  1  ^ 

Sqs.  (3-3-11)  end  (3-3-23)  provide  two  ordinary  linear  simultaneous 
equation*  in  the  unknown*  x  and  a. 

We  now  sake  a  transformation  in  the  independent  variable  by  putting 

* 

Z  =  *x  =  -e 


i 


(3-3-24) 
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then 

-  2  Jit  _  Hi  (3-3.25) 

»■»  ~  D*  ?  3* 

In  view  of  (3.3*24)  and  (3.3*25)  equations  (3. 3. 11)  and  (3» 3* 20) 
became: 

K.<L?  +  ,S4«'i  *-d.15  =0  (3-3.26) 

di  d*  cXl* 

4  2Ji  *<*15  +  4.^15  +4a  (3.3.27) 

oli1'  di 

We  first  examine  the  case  of  *■  ^  Z  .  For  this  range  of  values 
of  tA,  ,  (3.3.26)  and  (3*3*27)  yield  the  solution 

*i  -  ul  Ow)i  -(m+zH 

=  A(w)  €  *+•  6(*^)  •£  CW  €  -4-  El  f*)  €  (3.3.2o) 

_  ni  _  -Mi  (M-Z)l 

S  -  -  2A(*)  «  vi f3(**)e -lC(^)-e  (3.3.29) 

Transforming  equation  (3*3*2)  and  in  view  of  (3*3*10)  we  get 

K  6  =  *£  +  4*  -4-  I'E  (3*3.30) 

d * 
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hence  ^  is  found  to  be 


c  =  ^  Am*  Rm<  -  Ccm  «  *  EH *  (3.3.31) 


Two  of  the  unknown  constants  can  be  eliminated  from  the  condition 
that  *X  ,  and  hence  are  aero  at  the  inner 

boundary  i.e.  at  x  -  1  or  z  -  0.  Putting  *  «  0  in  (3.3.28)  and  (3»3*3l) 
we  get: 


A  +  ^+Ch-  ^  -^-O 


(3.3.32a) 


(3.3.32b) 


Expressing  C  and  E  in  terns  of  A  and  B  ve  get 


C  =-0--^A  + 


l/(- 


E  -  -  A  A 

1^ 


lA. 


(3.3.32c) 


(3"3*32d) 


Substituting  for  C  and  E  in  ( 3.3.S8)  and  (3«3*3l)  we  get 
.  (  Hi  ,  (w-i)i  -b+z'tl  ( 

r  :A»|  «  -  0-i)«  ~i  «  5 

a,  C-«  1  .  fa-1)*  ,  _fa-n)*l 

-£■«  J 


(3.3.33) 
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(u-iU  ,  -(«**■)*( 

e  ~  J 


(u-i* ,  .  -M* 

•e  -C+0« 


1 


(3-3-3^) 


Eq.  (3.3.33)  end  (3-3-3M  in  tense  of  the  independent  variable  x  become 


*  -  A  ^TLk_(l-£')-xM"-^  X  i  X -f^ i)x+lj(3.3.35) 

6  -A  *8{W--A  ST4*fji^>.3-» 


also 


S  ^ "‘V  irfi'-i-Z'i  (3-3.37) 

To  derive  the  displacements  we  note  eq.  (3.3.12)  and  (3. 3. 17) 
from  which: 


_ 

za 


(3.3.38) 


and 
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(3.3.39) 


Integration  of  eq.  (3*3.38)  will  Introduce  a  constant,  however, 
substitution  in  (3. 3. 12a)  for  ^  and  conparison  with  (3*3.37)  thews 
that  the  constant  wust  be  zero. 

Hence  on  integrating  (3*3*38)  we  find: 


II 

1 

-fc. 

S' 13 

'  K-»-l 

>  X 

l*-l 

X 

^  k-*-l 

k. 

p>  -  5 
6 

:  •-*- 
li- 

L  k—  * 

1 

k(u-0 

k(k-*-0  j 

u-i  , 

*.  _  J_  X 

\A- 


(3.3AO) 


Also  (3.3.39)  yields, 


-  =  Ajr.  x-  _i^ 


k-H  .  W-l 


(5*  ^  k (K+O  l*- 

S  2$-  'til 

G  i  k(k- 


~(u+-l 
X 


U.(k-H') 


1 


l-M*  H-l  —  (m-h)| 

\  — w  .  > 


X  +  - -  X  ^  —  x 

i')  k.(K-0  *A- 


( 3*3* **1) 


For  convenience  we  write  (3-3*  35)  i°  the  font: 

X  r  AM'P  «  +  ^  f  w 


( 3. 3.^2) 


Ik 


end  (3.3.36)  as 
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£  r  +  6w 

where 


Also  we  write  (3. 3.40)  as 


vl  3  _  4  K)  —  u  (so  _  U  (*) 

G,  Gi  2K 


and  (3-3.41)  *» 


IT 


-  AH  -  V  (*) 

G* 


-4  Sk>^-  V  M- 

&  2-* 


where 


•  X  I  K-*-l  ,  M— 1  I 

U  -  -L  X  __L^c  _  _  ^ 

>>A-  U  +  l  \A.  U.(U-M') 


( 3- 3*^3) 

(3.3.44a) 

(3.3.44b) 

(3-3.44c) 

(3-3-44d) 

(3.3.45) 

(3-3.46) 


(3.3.47*) 
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U  CO 

i*. 


I-U. 
X  -4- 


I 


k-l 


K-> 

x  _  x 

IA.(  U-l'') 


(3-3-47b) 


V  (*) 

m.  ' 


xM4l  i 

k(u+0  u- 


L-l 


■x 


u(u-*-0 


(3-3-47c) 


.  1— **•  .  m-i  -4M+r> 

V  M  =  -  x  +  — X  -t-—  x  (3-3-W) 

In.  u^(  K-0  *.(  U-l  )  l*- 


The  constants  A  CO  and  ^  (O  cannot  be  evaluated  since  nelgher 
the  stresses  nor  the  displacements  are  known  at  the  cylinder- shell 
Interface. 

For  the  case  of  *  «  1  equs.  (3-3-26)  and  (3-3-27)  become: 


ols,  _  d 


o  ,  oi*i  d't-t  _  <?  u  /V  _ « 

-  -4-  z  —  s  —  0 

d*  d?1  d* 


(3-3-48) 


CilS'  H-S,  =0 

ol%  *■  d?  d  i 


(3-3-49) 


where 


S  r 


s(?,o  %  =  *(*,0 


(3.3.50) 


Eq.  (3-3.46)  and  (3-3-49)  yield  the  solution  in  the  x  variable 


%  r  A  ,  x  -v  ^  1 


-3 


(3-3-51) 


Transforming  (3*3*2)  for  ■'t  ■  1  and  using  (3*3*iO)  we  get 


(3*3*53) 


and  substituting  for  S,  and  'Y ,  we  get 


-3  _  < 


<*,,  =  A  -x  *ft,  x+  E,-x  -fj-r. 


(3-3.5M 


Using  (3*3* 

38)  and  (3*3*39)  for  Vt  . 

1  we  get 

-  — 

2  Cx  l 

+c>} 

(3*3*55) 

O’ 

I 

?3A,xX-6*- 

E  (?«X*0-R> 

(3*3*56) 

r 

2G.  *■ 

*  J 

at  x  -  1,  'Y,  • 

o 

*• 

0 

o 

(3*3*57) 

hence: 

A  ,  4.  8(  -  O  (3*3.58) 

A  R  -vE_  _  p*  -r  -  o  .  . 

1  .  i  IT  •  (3*3*59) 


(3*3*60) 


1 1/~> 


(3.3.61) 


-  A.Cx-x’) 

(3-3-&) 

U  =  l!  £4  x-t-A^  -+*»<*.?©«%-.'''?  V  (3.3*63) 

'  ZG  L 

-R  J  (3.3.6U) 

1  2G  L  ’ 

We  finally  deal  with  the  case  of  n  »  0.  We  sake  use  of  (3*3.13) 


fe„-^0  =°  .  — °  + 

’  e  '  ol-f 

a°  -  0 

'V 

(3.3.65) 

=  i- 

(3.3.66) 

-v 

I  -  _  *1  I 

^  1  ’  e0 

A  a 

x 

(3.3.67) 

__  / 

50  -  -  2G  ^ 

(3.3.66) 

x 

w 


From  (3.3*1) 
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hence: 


c/  <*-»o 

2.  G 

(3.3.70) 

* 

O 

nr  0 

-  4;  - 

<£r  Ao 

“v” 

(3.3.71) 

1 

-  0 

(3.3.72) 

Let 


Then 


(3-3.73) 


-  A°  0  ^a.) 


(3.3.74) 


■ —  j  '-f  1  —  I 

AA0  —  Ao  ~~~ 

G» 


(3.3.75) 


Once  A( k )  and  B( )  have  been  found  from  the  stress  and  displace- 
»nt  continuity  conditions  at  the  cylinder- shell  interface  then 

t cx.osiu-xe  (3-3-76) 


HS  I 


<2  -  <2  +  21  <0  (*,  ^  Cos*© 


no  *» 

U  a  I 


(3.3.77) 
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(3.3-78) 

(3-3.79) 


r-  ^ 


The  equilibrium  relation*  for  a  abell  element  are: 


In  the  tangential  direction 


de 


(3.4.1) 
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Zb  the  noraal  direction 

d(2e  +  =  H  (o 

©  *  * 


and  the  sonant  equilibrium  ia 

Oe  r-L^s 

9  "*  de 


<3-k.  8) 


(3-4.3) 


Eliminating  Qe  from  (3.4.1),  (3.4.2)  and  (3*4.3)  we  get 
dropping  the  suffix  0 


dVI  ^ 

— .  -4-  — 

dG 1  dG 


+  „  a  _  I  cM 
dG  1  ^  ^  dG 


l  nj  -  -  k 
^  dea 


(3.4.4) 

(3-4.5) 

(3-4.6) 


The  force  displacement  relation  becomes 


N 

_  E..k 

j_  ^ 

(3.4.7) 

i-vr 

l  ^  dQ  J 

M 

0 

^  ol3i<s  d  & s  l 

(3-4.8) 

/»•* 

* 

1  7ex  ~~  di>  5 

where 


E.l3 

iz  O'V) 


P 


(3>-9) 
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Taking  Fourier  transform  of  (3.4 .4,5a,  6)  ve  get 


—  H_  Nj  -1-  H_I 

^  +  kJ  =  %  (0 

(3.4.4a) 

-  -  -  M 

(3.4.5a) 

-  M  -+-  M 

r  ^f3 

(3.4.6a) 

Also  taking  Fourier's  transform  of  (3-^*7)  and  (3.4.8)  we 

obtain 

NJ  =  C  [  k?4  *  ut  ^ 

(3.4.7a) 

M  =  ^P*[ 

(3.4.8a) 

where  p**.-  P 

C*-  E‘  ( 

’M 

s 

‘*V 

Substituting  for  |Vf  and  f\/  in  (3.4.5a)  we  get 
r  4WS  |  4  f  vC  n 


-*x\  -  j  (3.4.10) 

Similarly  substituting  for  M  and  Nl  in  (3.4.6a) 
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*vt jo  -  ^  P  ^  n1'  ^  ^  |  ^  ^  ^ 

or, 

^  -  ^*.4P*-*  C*^j  u^  ^u’P-^C*^  (3.4.11) 

For  continuity  of  stress  at  the  cylinder- shell  interface  ve  nust 
have  (Pig.  3-4.2): 


Figure  (3.4.2) 


jo(©)  z  -  <^(vt$)  -  WCe')-!-  uj  cjoiQ  (3.4.12) 

r  _  U)s  S<u€)  -  ^K,G)  (3.4.13) 

where  v#  is  the  veigbt  of  unit  length  of  shell  per  unit  length  of  arc. 


Bence: 
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-*(*»)-«*  ,  K  »  I 


-  I 


(3-4.14) 


(3-4.15) 


Also  continuity  of  displacements  at  ths  interface  requires  that 


U(^)=u4  u  f^)  =  as  (3.4.16) 

u’fO  =  0;  t>K)  =  (3.4.17) 


Case  of  <*-  £•  2  •  ^  )  ■  a. 

v  nr,  ' 

Fro*  (3-4.15),  (3*3.42)  and  (3.U.IO)  it  can  be  seen  that 


uJkC^'P*^  ^(cV)  (3.4.15) 


vv 

T? 


\  cj>  *  &  <j>  («-')  J-  =■  u^C*+iSD* |  ^•hs* 9  j  ( 3.4.19) 


Also  fro*  (3.3.45)  and  (3*3.46) 
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'll,  =  _  A  -  U  (*■)  _  8  !i  u  (<0 

5 


(3.4.20) 


if,  =  All  V  O')  -4-  6  —  \/  O') 
s  <S  "*■  &  *"• 


O.^.Sl) 


Solving  for  A  and  B  from  (3.4.20)  and  (3.4.21). 


A(' 


K.)  r  f  u  \/  7  -4-  J  (J/  l 

,y  |  1  2k.  ^  2  k  J 


(3-4.22) 


8 w  =  ^  fa  \/ "-+  s  u,L  | 

i,  [  5  s  J 


(3.4.23) 


where 


V  /  r  \/  (a)  f  U  (a)  \/  (*.)  -  U  («0  V  (<mI  (3-4.24) 
<k-  t  a*-  iw  ik  a**.  J 


/  C  /  —  | 

U  =  U  Co.)  }  (J  CM  V/  (ft.)  _  (J  C«0  V  (W )  (3.4.25) 

I*-  'K  [in.  '**  IK  tto  J 


etc. 


Substituting  for  A  and  B  in  (3«4.l8)  and  (3.4.19)  and  collecting 


tens  ve  get: 


U  f  &  I"  'V  (a*)\/(A_^l>  C^V^l 

J  L  2  K  IK  'll*.  2  K  J 


IT 

s 


(a-*)  ^  (*■)  L)  ia"(P  •+£*)|'*o(  3.4.26) 
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«.s  ^  ~  4? (<°  -  m-4,  c  *  | 


-4-  IT  {&  [<£ Jft.) (J,  V)  - 4?Jto  U/®*] - * P'“C  }=  ^  (3.4.27) 

From  equations  (3.4.26)  and  (3*4.27)  and  U"s  may  be  found. 

Substitution  of  (3.4.26)  and  (3-4.27)  in  (3-4.22)  and  (3-4.23) 
gives  A(^  )  and  B("-)  respectively. 


Case  of  n  ■  1 

Equations  (3.4.10)  and  (3.4.11)  become : 


\\  -  c*{  •+  P*{^,+U<A  (3-1,'26) 

_  ^  ^  I  ^  ^  / 

-  C  1  J-1-  P  )  (3-1*-29' 


hence: 

\  <3-'*-3o) 

The  relation  ^  satisfies  the  condition  of  vertical  equili- 

'1 

briua,  i.e. 


tt 


A 1  w  £)  d  (!) 


•  I 


cotO  d0 


-t, 


(3-4.31) 
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also  for  n  ■  1  fro*  (3-3-6)  and  (3-3-62) 


i  -  A,  fa-  i3}  (3-3-6l> 

g  =  (3-3.6S) 


In  view  of  these  and  also  (3-4.14)  and  (3.4.15) 


Q  r.  -  -A  ,  ^  3  )  -  (3-4-32) 

*» 

_  _A,  (a-<X,)  +  p-r|(a-£)-Mo  -  ^  (3-4-33) 


Equations  (3.3*63)  and  (3-3-64)  can  only  b«  consistent  In  view  of 
If 


W  -  2T’^1»0s-<-p  (-1,  -  )  (3-4-34) 

which  Is  true,  since  the  right  hand  side  represents  the  combined  weight 
of  shell  and  cylinder. 

Also  from  (3-3-63)  and  (3-3-64)  we  have 

JL  =  -  ^  ((X-O?)  -f  +  t> 

'  2&  l  ' 


(3-4-35) 
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\  A, a -  n  -  % ) 


Adding  tbeae  last  two  equations  we  get: 


Hence : 


A(  rfa^d2 


i 


also  fron( 3.4.32) 


Proa  (3.4.38)  and  (3-4.39), 


(3.4.36) 


(3.4.37) 


(3.4.38) 


(3.4.39) 


(3-4.40) 


Otherwise  by  eliminating  (  utS(  4  CTSi  )  we  can  determine  ^  i.e. 
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or 


-4,(a-d3)  ^s  +  i;  7^  (3 AM) 


Naturally,  for  the  nosant,  U.^(  -  cannot  be  determined  since 

it  represents  a  rigid  body  displacement.  Finally,  hcwerer,  US)  will 
be  dete rained  from  the  condition  that  U(  >  0  at  0*0. 


Case  of  n  ■  0 


Froa  (3.4.11) 


=  a5.C 


(3>.W) 


also  in  view  of  (3-4.14) 


b  = -4  _ 

-n"! 


(3.*.*3) 


b  -  _  A(\-±  — 

to  0  a1 


(3.^. 43a) 


Hence: 


/•  * 

vtIoC  = 


(3.4.44) 


But  froa  (3-3.76) 


.  A  li.- 

0  6i  o- 


(3.3.75) 
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or 


Ao  = 

Substituting  for  Aq  in  (3.4.44)  we  get 


W 

IT 


u$o  =  _  jV+G  <*-''>]' 


also  fron  (3.4.45) 


or 


A  aNA/k  fc*-*-  6(a2-0  \ 

0  '  TT-r,  C  * 


(3^.45) 


(3  AM) 


(3.4.4?) 


This  complete*  the  formal  solution  to  the  elastic  problem. 


3.?  The  Associate  Viscoelastic  Problem 

To  obtain  the  solution  to  the  viscoelastic  problem  ve  make  use  of 
the  Integral  form  of  the  constitutive  relation,  l.e. 


G  (4  d% 

3* 


(3.5.1) 


-«0 


Upon  taking  Laplace  transform  with  respect  to  time  (2.3.I)  becomes 


(3-5.2) 


10? 


where  the  dash  designates  a  transformed  quantity.  In  riew  of  (3.5*2) 
It  is  evident  that  the  solution  to  the  viscoelastic  problem  can  be 
obtained  on  substitution  of  G  b  for  G  in  the  relevant  equations 
and  taking  the  inverse  transfora  C  l  3  • 

The  numerical  analysis  of  the  problem  consists  in  calculating 
the  displacements  of  the  cylinder  both  at  the  inner  and  outer  radius. 
These  results  will  be  presented  at  a  later  date. 
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APFCTDIX  I 


Hu— rlcal  Applications 

Infinite  Viscoelastic  Slab 

To  test  their  accuracy,  the  — thods  outlined  in  Chapter  I 
were  applied  to  the  problea  of  a  viscoelastic  slab,  which  has  been 
solved  exactly,  within  the  accuracy  of  the  mi— rical  computations 
employed,  by  Mukl  and  Sternberg 

The  slab  is  of  infinite  extent  and  finite  thickness  2a,  bounded 
by  the  planes  Z  -  -  &  •  Initially  It  is  at  unifora  temperature, 
8o°C,  when  suddenly  the  temperature  of  the  faces  changes  to  110°C  and 
re— ins  constant  at  that  value. 

The  solution  of  the  problea  reduces  to 

-  Si  =  =  °  -  *  -  f<V  H 


y 


-  e 


71 


-  £ 


~o 


= 


-Gy  =  O 


(A. 1.2) 


where 


is  found  froa  the  integral  equation 


J_  (  £(£-£'>?**  d*  (A.1.3) 

3<\  S  ^  l  Z* 

'  O 


In  Ref.  [V  J  the  problea  was  solved  by  referring  (A.1.3)  bo  the 
plane  and  taking  Laplace  transform.  However,  in  view  of  the 
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realistic  nature  of  E(t),  finding  the  analytic  representation  of  the 
latter,  taking  Laplace  transform  of  (A. 1.3)  and  Inverting  was  a  very 
tedious  procedure. 

Here  we  obtain  a  solution  at  the  middle  plane  of  the  elab  first 
by  reducing  (A.1.3)  to  a  set  of  linear  algebraic  equations  by  the 
method  described  in  Chapter  I. 

Let  be  piecewise  linear  in  t  i.e.  in  the  interval 


Also  let 

4 ra^  (A,1,5) 

Then  in  matrix  form  (A.1.3)  becomes 

WH  L£]  (“"l  (A-1-6) 

where  F  is  a  known  function. 

From  (A. 1.6) 

M  =  C '  { f}  (A.1.7) 

where  [£]  "  LEJ  (A.1.8) 

■*  =LHJta}  (A,1,9) 


Referring  to  the  data  of  Ref. 


for  the  Material  properties  we 


can  find: 


2.24 

2.01 

2.24 

1.79 

1.88 

2.24 

1.52 

1.61 

1.69 

1.29 

1.33 

1.39 

1.00 

1.04 

1.08 

.80 

.80 

.80 

2.24 

1.44 

2.24 

1.11 

1.15 

2.24 

.80 

.80 

.80 

2.24 


.10 


10 


(A. 1.9a) 


R 


1.2933 

1.2632  1.2933 


1.2344 

1.2462 

1.2933 

1.1991 

1.2019 

1.2213 

1.2933 

I.1689 

1.1742 

1.1820 

1.1886 

1.2933 

1.1397 

1.1362 

1.1414 

1.1454 

1.1506 

1.2933 

1.1048 

1.1048 

1.1048 

1.1048 

1.1048 

1.1048 

.10 


10 


(A.1.9b) 


1.2933 


(A.  1.9)  being  the  solution  to  the  problea.  The  results  obtained  by  this 
■ethod  have  been  plotted  for  comparison  against  the  "exact"  solution 
obtained  in  Ref.  [jj  -  See  Pig.  (A.  1.4)-  and  they  compare  very  favorably, 
in  view  of  the  fact  that  only  a  desk  calculator  was  used  for  the  numer¬ 
ical  computations. 


Ill 


The  iterative  aethod  vu  alto  used  to  solve  (A.  1.3)  again  at  the 
middle  plane  £  -  O  .  The  zero ' th  approximation  to  the  solution  was 
found  by  putting  K  =  o o  .  Then 


t 


=  -2*. 


© 


Eff  -§')  d*i 


(A.1.10) 


Iteration  was  then  carried  out  on  the  basis  of  the  foxmila 


i 


E(f d*  r 


i 

z  -  z«.f  EXf-fO^d* 
)  a* 

o 


(A.l.ll) 


were  in  iterations  subsequent  to  the  zero' th  the  actual  value  of 
was  used.  In  Fig.  (A. 1.4),  it  is  shown  that  after  two  iterations  the 
solution  obtained  was  practically  -  within  error  of  numerical  computa¬ 
tions  Involved  -  identical  with  the  "exact  solution"  cited  above. 

The  above  iteration  procedure  is  in  fact  convergent  as  can  be 
seen  by  the  following  simple  proof. 

Writ.  E(f-f')  i  (*.1.12) 


where 


£(f-V>  jf  i 


(*.1.13) 
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O 

A 

0° 

h 

v-£- 

V 

n 

? 

ot 

1 

II 

uis 

(A.1.14) 

Equation  (A. 1.6)  becomes: 

t 

<* 

4-  ^  r  scf-?')  d*  -  ~  i-  (t-) 

■’o  &  ^ 

where 

i 

F(tr)  =  fl*.  d* 

J0  a't 

(A.1.15) 

Furthermore  F£c)  la  bounded  In  every  Interval,  o^'t  $  t  f 

See 

Fig. 

(A.l.U). 

Hence 

|  F(t)  1  £  M 

(A. 1.16) 

Then 

(A. 1.17) 

t 

=  -  FCO  -  ^Ecf-r}^  d* 

(A.1.18) 

and 

1^,1  * 

(A. 1.1Q) 

or 

Kx,  )  ^  M  +  jM 

(A. 1.20) 

I  *  M4JM4 


also 


(A. 1.21) 
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and 

£  M  +'f)  (A-1-28> 

.  k. 

The  series  - i-}  is  convergent  for  I  and  equals  — 

Hence 

|£x  I  <  J±  (A.1.23) 

'-5 

and  the  iteration  is  convergent.  Also 

?  Mr  —  ° 

k.  — #»  oo 

Hence  solution  can  be  approached  as  closely  as  ve  please,  by  a 
sufficient  number  of  iterations . 


Figure  (A. 1.1)  TIME  (hour*) 

TER3I0H  MODULUS  E  FOR  METHYL  METHACRYLATE 


2 


TIME  (hours) 


3 


Figure  (A. 1.2) 

AT  MIDDLE  FLAKE  OF  SLAB 


lIjn 


^  (hours) 
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TIME  (hour*) 

Figure  (A.1.4) 

STRESS  HI  MIDDLE  FLAKE  OF  THE  SLAB 


11 B 


Viscoelastic  Sphere  With  a  Step  Rise  la  Surface  Temperature 

The  ine  methods,  under  sore  unfavorable  circumstances  of  non- aero 
initial  conditions  and  rapid  changes  in  teaperature,  are  now  used  to 
obtain  the  hoop  stress  at  the  surface  of  a  solid  viscoelastic  sphere , 
subjected  to  a  step  teaperature  rise  at  the  surface,  it  is  found  that 
in  applying  the  first  asthod  soae  care  aust  be  taken  near  t  »  o.  Other¬ 
wise  both  aethods  give  results  which  are  in  excellent  agreement  with 
the  exact  solution  of  Ref  .C  73  . 

Mukl  and  Sternberg  £73  obtained  an  exact  solution  for  the 
stresses  by  assuming  that  the  material  has  an  elastic  bulk  modulus. 

Using  their  notation,  can  be  obtained  by  solving  the  follow¬ 

ing  integral  equation, 

t 

dY  (A.1.25) 


whilst  on  the  surface  of  the  sphere 


£  K ,  t)  Z 

9 


_  1  S  (•*,  .  t) 

2  ’ 


(A. 1.26) 


When  (A.1.25)  is  referred  to  the  (  -r  f  )  plane  it  becomes 


f  4  ^  A 

swt)+l  (©MVl-d?'  c  oie'  (A.1.27) 

’  '  1*  S.  35'  Jf' 


where 

„  A 


^  I  * 


(A.1.26) 


Nov  taking  Laplace  transform  of  (A. 1.26)  with  respect  to  ^ 
,  one  obtains: 


or 


* 

S 

or 


1  <  S  [?  k 

3^+26  (=» 


vbere  f?  is  then  obtained  by  Laplace  inversion. 


It  is  shewn  however  in  the  General  Introduction 


(A.1.29) 


(A. 1.30) 


G  -  ^  K  ^  (A.  1.31) 

9K-f»£ 

where  E(t)  the  tension  Modulus  of  the  Material. 

On  substitution  of  (A. 1.31)  in  (A.l.30)  we  obtain: 


V".  P 


2  <  £  ,  f  .  , v 

- =r 


(A.l.32) 


On  taking  inverse  Laplace  transfora  of  (A. 1.32)  and  reverting  to 

the  (  *,  t"  )  plane  (A. 1.32)  becoass: 

t 
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On  tbs  other  hand  fro*  (A.1.30): 

t 


Skt)  -  f 

-  '  )  2*1 


(A.1.34) 


naturally,  since  R  is  known  a  solution  for  can  be  obtained 
fro*  (A. 1.34)  by  staple  quadratures.  Such  a  solution  provides  a 
basis  for  c caparison  with  the  two  suggested  Methods  of  solution  of 
(A.1.33) 

Solution  of  (A.1.33)  by  reduction  to  a  set  of  algebraic  equations. 

We  proceed  to  solve  (A.1.33)  on  the  surface  of  the  sphere. 

In  Chaf>.r  it  was  suggested  that  equations  similar  to  (A.1.33) 
could  be  solved  by  assuming  a  piecewise  linear  variation  of  the 
unknown  function  with  tine,  assuming  that  the  function  was  zero  at 
"t  -cm-  .  nils  is  not  the  case  here,  however,  by  virtue  of  the  right 
hand  side  of  (A.1.33)  which  is  different  froa  zero  at  1  =  o-+-  .  On 
the  other  hand  it  is  convenient  to  solve  an  auxiliary  equation  in 
which  the  right  hand  side  is  ^  £  .  i.e. 


^(ew-L  f 


(A.1.35) 


r  VC.  . 
>  =  IO  C 


On  the  surface  of  the  sphere  we  have  fro*  the  example  of  Ref. 
,  and 


E -If'’.,*}]  I  £[loVSt-t')j  =  £(1-0 


(A.  1.36) 
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If  we  look  at  (A.1.35)  la  the  (  f  )  piano,  then,  after  Laplane 
tranafora  and  in  view  of  (A.1.35)>  and  (A.  1.33)  v«  get 

i 

s  A  jAt  (A.1.37) 

>3  d* 

o 

after  comparison  of  (A.1.37)  with  (A.l.30)  it  transpires  that 

=  j?  [  (A.l.38) 

Since  R  is  knovn,  sn  estlaate  of  the  accuracy  of  the  aethod  can 
he  aade  on  the  basis  of  solution  of  Eq.  (A.1.35).  Again,  of  course, 

^  o 

Hence  following  Ref.  Q12]  : 

r  H  (t)  ^(0  +  )  (t)  (A.1.39) 

where  H(t)  is  the  unit  step  function,  and 

Substituting  (A.1.39)  in  (A.1.35)  and  in  view  of  (A.l.36): 

'vlVo  H  (e^  +  \P*n  +  EL(tUi  (e.(+-t )*-d*  =1  £  (to  (a.i.4o) 

T  7  3*  W)  2*  3 

9* 

Putting  "t  =  ©+  we  find 

^(.+)  r 

In  view  of  (A.l.Ul),  (A.1.40)  becoeas 


( A.l.Ul) 
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■  -[  Eft-t^a*  '  \  (a.i. 

3*1  3t  E.(e>+)  * 


Eq.  (A.  1.39)  1*  »  necessary  step  and  the  emission  of  the  initial 
conditions  would  introduce  significant  error  in  the  solution  of  (A. 1.35) • 
The  right  hand  side  of  (A. 1.42)  is  now  a  known  function  vhlch  we 
denote  by  g(t)  i.e. 

Aw 

«(n  =  -  ^(ovi-T,  _  1  (*.x.*3) 

J  I 

and  (A. 1.42)  becoaes 

#  *  * 

^V)  -4.  J  ECfr-O3-^  di^t  =  <jte  (A.  1.44) 

Because  of  the  nature  of  E.U-)  Fig.  (A. 1.8)  g(t)  has  a  very  high 
gradient  at  t  ■  o. 


Figure  (A.1.5) 
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This  is  a  potential  source  of  error  for  a  finite  initial  interval, 
unless  the  first  intervals  are  taken  extremely  snail.  To  avoid  such 
a  lengthly  numerical  computation  we  approximate  here  the  high  gradient 
of  g(t)  at  t  «  o+  by  a  finite  step  at  t  -  o  as  shown  in  Fig.  (A.1.5). 

Thus  again  we  write 

<^{t)  -  +  H  O-'*  +  (A.1.U5) 


and  consequently  we  can  write: 


¥■  v 
^  (£) 


(A. 1.46) 


Substituting  (A. 1.46)  and  (A. 1.44)  we  obtain 


and 


* 

vf>  (o- M-) 


*3(0+-^ 

|  4-  E.Co+1 

Tk 


in 


(< 


(A.1.47) 


(A. 1.48) 


The  right  hand  side  of  (A. 1.48)  is  now  a  well  behaved  function 
and  (A. 1.48)  can  be  solved  by  piecewise  linearization  of  over 

snail  finite  intervals. 

Naturally: 

*  / 

-^(0+)  +Npco+-o  +  it) 


(A.1.49) 
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Solution  of  ^(f-)  baa  been  obtained  over  a  total  tine  of  .20  hours, 
in  the  following  intervals: 

o,  o+ ,  o++,  .01,  .02,  .04,  .08,  .12,  .20 

^ ,  which  is  Identical  toR  ££(*,,{-  )J  ,  thus  obtained  is  ccsipared 
against  the  exact  solution  of  Ref.[7 J  •  The  agreement  between  the  two 
solutions,  shown  in  Pig.  (A.  1.10)  is  exceedingly  good. 

The  hoop  stress  was  then  calculated  frae  (A.  1.37)  and 

(A.  1.26)  and  is  coopered  with  exact  from  (A. 1.34)  in  Pig.  (A.l.ll) 

Again  agreement  is  very  satisfactory. 

Solution  of  (A. 1.33)  by  an  iteration  technique. 

Since  ’'■f'Vo-O  has  been  found^eq.  (A.1.41),  we  solve  (A.1.44) 

Instead.  In  other  words  in  the  xeroth  approximation  ^  is  considered 
infinite  for  all  values  of  time  except  at  t  »  o+.  This  apparently 
artificial  physical  consequence  can  be  appreciated  if  one  looks  at 
N^(t). 

Prom  (1.5*0 

VCO  = 

Hote  that  \?  (O 


z. 


E.Ct-’t 


(a.1.50) 


— *■  as  H  -*■  °0  ,  See  Fig.  (A.1.6) 
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i 


Figure  (A. 1.6) 

Hie  above  values  of  K  imply  that  at  t  ■  o  but  \?  -  y 

for  t  y  O  ,  In  other  vorda  ve  have  approximated  \?(t)  by  the  thick 
line  shown  In  Fig.  (A.1.6). 

Is  then  found  from  the  recurrence  relation 
*  + 

^  -  *\(t)  (A.  1.51) 

^  2*  * 

for  n  ■  1,  2  ...  where 
•if- 

^  (A.1.52) 

O 

Values  of  obtained  by  this  method  are  also  shown  in  Fig. 

(A.  1.10),  and  these  too  compare  favorably  with  the  exact  solution.  Hie 
hoop  stress  is  also  found  and  plotted  in  Fig.  (A.  1.11),  with 

good  agreement. 
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Stresses  In  Viscoelastic  Solids  du«  to  Cyclic  Teaperature  Histories 

The  effect*  of  cyclic  temperatures  on  the  stresses  In  viscoelastic 
bodies  ere  not  veil  understood,  and  this  Is  not  surprising  In  view  of 
the  analytical  difficulties  Inherent  In  the  Investigation  of  this 
p rob lea.  The  critical  effect  of  teaperature  on  Material  properties, 
renders  the  p rob lea  non-linear  in  so  far  as  the  principle  of  super¬ 
position  no  longer  holds  with  regard  to  separate  teaperature  histories. 

.0^  fA\  .  —  «  . 

That  is,  If  T(-x  produces  £  .  ("x  t- )  and  \  ,  t- > 

produces  (x  J-')  then: 

(A'  *  *  —(ftt  . 

T  (%,(•')  +  T 

does  not  give  rise  to 

(A>  ^ 

The  superposition  principle  will  hold  only  in  the  particular 

_ (a-)  (a) 

case  where  \  ^  T  are  siaultaneously  taken  to  affect  to  Material 

properties  but  their  effects  on  the  strains  are  considered  separately. 
In  such  a  context,  however,  superposition  will  be  of  little  value. 

To  investigate  the  phencnenon  at  all,  we  take  the  simple  case  of 
the  viscoelastic  slab,  and  subject  its  surface  to  cyclic  teaperature 
changes  shown  in  Fig.  (A. 1.7). 


We  •xniM  the  effect  of  the  two  cycles  shewn,  on  the  street  »t 
the  surface  of  the  Infinite  slab,  when  (h^  ■  50°C,  and  the  datum 

temperature  is  70°C. 

Moreover  we  assume  that  the  slab  is  aade  of  LPC-5^3A  propellant 
for  which  the  shift  factor  is  a  well  town  function  of  temperature, 


see  Ref.  ,  p.  55* 

For  this  particular  material  it  is  found  that  the  shift  factor 
is  75  when  (3  *  50°C  ,  and  when  (3  “  -50°C  or 


f-{ 


1st  ©=  5o*  C 
L  t  &  -  -  yoc 


(A.1.53) 


IS 


In  the  notation  of  Ref.  ,  on  the  surface  of  the  slab 

I 

o  r 

Fig’s  (A.  1.12)  and  (A.1.13)  ahcw  ^  as  a  function  of  t  and  (© 
as  a  function  of  ^  .It  is  easily  seen  that  ©  can  be  expressed 
in  the  following  form. 


©  =  +  J  (A. 1.55) 

and  consequently 
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where  it  the  delta  function  at  £  -  £  .  Substitute 

(A.l.56)  is  (A. 1.5**)  end  integrating  ve  get 

r  <A-l,57) 

From  (A.  1.57)  ve  can  easily  find  as  a  function  of  time  for 

two  temperature  cycles,  see  Figs.  (A.  1.1**)  and  (A. 1.15). 

It  can  be  seen  that  in  both  cases  increases  asymptotically 

to  values  which  are  almost  twice  the  maximum  for  identical 
temperature  cycles  applied  to  an  elastic  slab.  This  phenomenon  will 
take  place  irrespective  of  the  period  of  the  cyclic  variation,  and  is 
due  to  the  difference  in  the  relaxation  rates  at  the  lower  and  higher 
temperatures . 

Though  this  result  is  important  it  cannot  be  readily  generalised 
either  for  more  complex  temperature  histories  or  for  more  intricate 
boundary  value  problems.  It  does  however  bring  out  a  phenomenon  that 
must  be  considered  in  the  design  of  solid  propellant  configurations. 


129* 


Flgux*  (A. 1.8) 


TIME  (HRS) 
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Figure  (A. 1.9) 
FUNCTION  h(t)  n.  TINE 
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SPHERE  TIME  DEPENDENCE  07  <p##  AT  SURFACE 
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Figure  (A. 1.12) 
Cycle  (*) 


Figure  (A. 1.13) 
Cycle  (b) 
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APPENDIX  II 


Error  Analysis  of  Approximate  Solution!  Developed  in  the  First 
Chapter 

Introduction 

In  Chapter  I  it  was  shewn  that  the  solution  of  viscoelastic 
boundary  value  problems  associated  with  the mo rheologic ally  simple 
viscoelastic  solids  in  non-uniform  transient  temperature  fields 
reduces  to  the  solution  of  a  Volterra  integral  equation  of  the  second 
kind.  In  some  cases  [7]  the  integral  equation  can  be  expressed  in  a 
convolution  fora  in  the  reduced  variable  and  can  then  be  solved  by 
taking  Laplace  transform. 

Where  this,  however,  is  impossible  or  where  the  relaxation  moduli 
of  the  material  are  in  the  form  of  experimental  curves,  which  is 
invariably  the  case,  Volterra  integral  equations  can  be  solved 
approximately  by  reduction  to  a  set  of  algebraic  equations  as  vas 
shewn  in  Chapter  I.  The  solution  of  these  equations  is  easy  because 
the  matrix  of  the  coefficients  is  triangular. 

In  this  report  we  investigate  the  error  inherent  in  this  procedure 
by  finding  first  an  auxiliary  solution.  Upper  and  lewer  bounds  to  this 
solution  are  then  established,  and  these  are  utilised  to  establish  an 
estimate  to  the  aforementioned  error. 

Upper  surd  Lower  Bounds  to  Solutions  of  Volterra  Integral  Equations 


Consider  the  integral  equation 
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v 

(j>(tf  +  3|  I"  =[(fc) 


(A. 2.1) 


where 


b 


d  x 


(A.2.2) 


and  ^  la  a  positive  real  nuaber,  and  KW  la  an  lntegrable  non 
negative  function. 

The  symbol  ^  denotes  the  reduced  tine  variable  given  by  the 
relation 


f 


-t 


T  a(*,t-)dt 

J  o 


(A.2.3) 


and 

{{H*, H) 

CUX.t)  -  €  (A.2.4) 

where  (T^  is  the  shift  function. 

Whatever  the  fora  of  T(x  ,  H  f  tor  any  "X  ,  a(x  ,  O  is  a 
non- negative  function  of  ~t  and  hence  ,  t"')  is  a  aonotonically 
increasing  function  of  ”t 


Eq.  (A.2.2)  aay  be  written  in  the  for* 
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K*(t  <fc')  -  £  Q  ^  t  )  ~^v«'  A  (a. 2. 


5) 


^ ^  and  assuming  appropriate  values. 

Examining  a  typical  tens  of  (A. 2. 5)  we  deduce  the  following 
properties  for 


K«\t)  -  k; 


(A.2.6) 


For  any  "t  ,  K'Ct,^  is  a  aonotonically  increasing  function 
of  At  ,  since  G  (f)  is  a  aonotonically  decreasing  function  of  t  , 
this  being  the  property  of  the  relaxation  modulus  of  linear  viscoelastic 
solids .  Also 


KVt,  o)  -  K» 


(A.2.7) 


where  K*  ^  ) 


is  aonotonically  decreasing  in 


-t 


Fig.  (A. 2.1 


t 

Figure  (A. 2.1) 


•t 


Also 


fco  = 


1UO 


(A. 2. 8) 


where  H  is  the  Heaviside  unit  function.  Without  loss  of  gener¬ 
ality  we  also  assuae  that  -^C°)  ,  hence  C°)  = 0  ,  because 

the  case  of  -j-0 0  ^  o  can  be  reduced  to  the  previous  case  as  will 
be  shown  later. 

Eq.  (l)  can  be  reduced  to  the  standard  fora  of  the  Volterra  inte¬ 
gral  equation  of  the  second  kind  by  integration  by  parts,  i.e. 

t  * 

3 ^  (A.2.9) 

O  O 


or 


cp(t-)-  ^  K  (t.OcjWci* 


(A.2.10) 


where 


Kfr.O 


9 

9* 


(A.2.11) 


3 


/ 


3 

I  +X  'Co 


i +x 


(A.2.12) 


Froa  (4.2.5) 
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(A.2.13) 


Now  since  Ckl't')  is  always  positive  and  G*(t)  is  aonotonically 
°  & 
decreasing,  Q  is  non-positive  and  hence:  is  non- negative. 

Le—a.  Subject  to  the  above  restrictions  on  *<(*,*)  ,  4)(0  is 
non-negative  and  non-decreasing  if  |((")  is  non-negative  sold  non-decreas¬ 
ing. 


Proof: 

Eq.  (A. 2. 10)  is  the  standard  fora  of  the  Volterra  Integral 
equation  of  the  second  kind. 

From  the  theory  of  integral  equations  it  is  proved  [  ll) that 


where 


* 


K  Ct.O 


is  defined  by  the  recurrence  relation 


*  f  *  * 


(A. 2.14) 


(A.2.15) 


and 


*  * 

K(  -  K(t.O 


(A.2.16) 
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*  * 
la  known  as  the  rnth  iterated  Kernel.  Since 


m. 


is  non-negative  all  the  iterated  Kernels  are  non- decreasing  in  *t 
Also  since  ^  /  is  a  positive  Member  it  follows  that  if  £(f)  is  a 
non-decreasing  function  dj)  (0  is  also  non- dec  re  as  ting. 

We  shall  utilize  this  important  fact  in  what  follows.  In  general 
will  be  an  arbitrary  function  so  we  seek  the  solution  for  (£>(♦") 
when  is  a  simple  non-decreasing  function  such  as  H  fr)  and 

use  the  solution  to  construct  another  solution  for  arbitrary  (V)  • 

Reformulation  of  the  Problem 

Given  (A. 2.1)  where  -j~(F)  is  an  arbitrary  function  we  seek  the 
solution  to  the  equation 
+ 

\p(H  +  }  [  -  H(tO  (A.2.17) 

Jo  ^ 


then 

<J)^)=^  (A.2.18) 

O 


Figure  (A. 2.2) 


Proof: 


{(0  -  Af(  H (*-+,)  +  £ H(+ '<*)-•-  ••••-<-  A f  (A.2.19*) 


See  Pig.  (A. 2. 2) . 

In  view  of  (A. 2.18) 


d)(0  =  Af  4  * ■■  +  A^MJ(A-2-19'>) 


In  the  limit 


4>(0  = 


d'c 


(A. 2. 19c) 


However 


O  bo  we  reduce  the  solution  to  the  standard 


case  by  writing 


Ap  ((-)  =  l-j  (t)  -+  A^  (f-) 


(A. 2. 20) 


where 


4*  (o)  =  o 

Then  substituting  for  4(0  in  (A.2.18)  we  get 


T  f 

+  (A.2.21) 


Putting  we  get 


+  -  1 


(A.2.22) 


hence 


^  - 


(A.2.23) 


i+><= 

and  substituting  for  'nI'  in  (A. 2. 21)  we  find 

1  o 

-t 

V(e>  +  >  f  K(+,  =  3_)V0-krfn} 


(A.2.24) 


and  the  right  hand  side  ie  now  a  aonotonically  increasing  function  of 

t  . 


Let 


•  +  ^ 


(0 


(A.2.25) 


then  (A.2.21*)  becomes 


-4 


-jr 

j  d't  =  <j(^ 


(A.2.26) 


and  'Vi/  is  a  aonotonically  increasing  function. 

Solution  of  (A.2.26)  by  Reduction  to  an  Algebraic  3et  of  Equations 
We  divide  the  range  of  integration  into  sasdJ.  intervals: 


1*5 


.  .O 

In  the  first  interval  ve  get 


By  the  mean  value  theorem 

<< 


where 


O  V<  »)  <  -t 


Letting 


V<*;)  H  ^ 


(A.2.27)  becomes 


S;  -*  -  3 


In  view  of  the  properties  of 


Kc+.-t) 


(A. 2. 31)  nay  be 


(A.2.27) 


(A. 2.26) 


(A.2.29) 


(A.2.30) 


(A.2.31) 

written  in 


terms  of  the  inequalities 
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>,  9, 

hence: 


^  <:  ^  $ 
I  +}  K® 


j>. 

l  +  ^.o 


In  the  second  interval  ve  get 

t  -3. 


<  -t 


1 


We  also  have  the  inequalities 

£  K'fM,')  *  <|4-.V> 

<(+„+,■>  ■?  ?  <■ 

*  .  *  . 

^  ^  >  o 

2  » 


(A.2.32) 

(A. 2. 33) 

(A.2.34) 

(A.2.35) 

(A.2.36) 

(A. 2. 37) 

(A.2.38) 


(A.2.39) 
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In  view  of  inequalities  (A. 2. 37),  (A. 2. 38)  and  (A.2.39),  (A.2.35) 

becomes : 


-v-  .  * 


(A.2.40) 


t*  +  3  <u,  ^ ^ 


(A.2.41) 


or 


><> 


^  „l)  £  "  '  (*.2.42) 


l  ^ 


z> 


and  aimilarly  for  other  intervale. 
It  follows  that  if  we  choose 


•  .  t- 


(A.2.43) 


we  get  a  lower  bound  to  ^  which  we  denote  by  and  by 

choosing 


V  -  fi- 


(A.2.44) 


which  we  denote  by 


we  get  an  upper  bound  to 

Hence  /\|^  is  bounded  by  the  solutions  and  ± 


or 


(A.2.45) 


Then  since 


/4'G-fc'> 


is  monotomic ally  decreasing  in  't 


making  use  of  Bonnet's  second  mean  value  theorem 


1  T  i k«x 


where 


Hence  if  ^ 


some 


1, 


o  <  t 

denotes  the  error  in  the  solution  of 


<  1;  ?  t; 


then 


z 


llAA> 
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(A. 2. 46) 


(A.2.47) 


and 


(A. 2. 48) 


for 


(A.2.49) 


5  'i-i)L.Jf(V)Lv 


(A.2.50) 
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where  again  o  ^  ^  <  "t 

Matrix  Formulation  of  the  Problem 

Coneider  (Fig.  A. 2. 3)  some  which  is  a  non- decreasing 

function  of  "t 


Figure  (A. 2. 3) 

ie  capable  of  two  approximate  step-wise  representations. 
One  is  the  full  line  and  the  other  the  dotted  line,i.e. 


or 


7  Hit-wii,., 


(A.2.51) 


^  --  NU-C.Ul, 

"vr  o 


(A.2.52) 


Substitutions  in  (A. 2. 27)  show  that  in  fact  (A. 2. 51)  yields  the 
upper  bound  solution  for  ^  and  (a.2.52)  the  lower  bound.  Thus  we 
get 
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J 


I  “i*  1; 


or  in  matrix  form, 

[H]f  + 


where 


(  i  I 


and 


'lO  ’'ll 


^?o  ^ %» 


< 

*«,»- 1 


From  (A. 2. 5^) 


fAfJ  {9} 


(A.2.53) 


(A.2.51*) 


(A.2.55) 


(A.2.56) 


(A.2.57) 
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and 

hit  -  [«][[H]*>N]'h! 


(A.2.58) 


Similarly 

W,j  -[h][lhJ^LM]  fi) 


(A.2.59) 


where 


CkJ 


Kl( 


k',* 


te"'  h  t 


(A.2.60) 


Naturally  a  linear  variation  in  ^  as  indicated  in  Ref.  1 
will  yield  a  solution  lying  between  the  two  bounds. 

Example:  (From  Appendix  i) 

Consider  the  slab  problem  where  the  relevant  Integral  equation 
to  be  solved  is 

i  * 

dx  +j_  f  e (?-?'> 

3K3.  d't  J0  d't 
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In  this  case 


,  and 


(A. 2. 62) 


We  first  solve  the  auxiliary  equation 


Aj;(f)  +  _L  f 


(A.2.63) 


and  obtain  upper  and  lover  bounds  to  .  by  first  finding 

as  Indicated  in  the  previous  section. 

It  Is  inmediately  obvious  from  (A. 2. 63)  that 


^  (pA  ^ 


£.  (•>+) 
=JK 


(A. 2. 64) 


A  solution  for  is  obtained  in  the  range  °  -2>  $  '7  hours 

and  the  range  of  integration  is  divided  into  the  following  intervals 
(-00,0+)  (o+,.l)  (.1,24)  (.24, .4)  (.4,-5)  (.5, .6)  (.6, .7)  hours. 

Under  the  above  conditions  ve  obtain  the  following  values  for  /'^ 


.01801 

•03591 

.05863 

.07942 

.10629 

.12822 


.01845 

.03676 

.06031 

.08160 

.10985 

•13197 


•77321 

.79122 

.80913 

.83184 

.85263 

•87950 

.91043 


.77321 

.79166 

•80997 

.83352 

.85481 

.88306 

.90518 
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For  numerical  details  of  the  analysis  see  12^  • 

Froa  the  table  above  and  from  Fig.  (A. 2. 5)  it  can  be  seen  that  /'V 


lies  within  close  bounds  the  maximum  difference  between 


and 


being  less  than  0 


Fig.  (A.2.6)  shows  the  exact  solution  lying  between  the  upper  and 
lower  bounds.  This  is  because  of  the  form  of  ^(t)  in  (A. 2.1).  In 
general  however  this  may  not  be  the  case. 

On  the  other  hand  ^  helps  to  bracket  the  error  and  guides  the 


choice  of  the  size  of  the  intervals  of  integration. 
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HOURS 


Figure  (A. 2. 5) 


fT.-T.)E 
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APPENDIX  III 


Proof  of  Convergence  of  the  Iteration  Solution  Obtained  In 
Chapter  I. 

We  finally  present  a  proof  of  the  convergence  of  the  Iteration 
solution  to  the  cylinder  problem,  formulated  in  Chapter  I.  Die 
solution  is  obtained  by  perturbation  about  the  uncompressed  state. 
Let  us  expand  these  stresses  and  strains  into  a  series  of  l/K: 


£ 

9 

£ 

'y 

(z 

6 

-y 

6 


oC 


(A. 3.1) 


r  Ca  4-  &  -4-  ^ 

&  i 


where 


(A. 3. la) 
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Since 


€ 


we  obtain 


*.  +  l 


£o 

=5* 


«o 


Comparing  power*  of  l/K: 


eo  =  3*0<s> 


In  general 


e  r  3^o0 


Substituting  for  t  in  terns  of  ^ 


fc  =  —  -  (Vt9)  - 


&u-' 


3^ 


K 


(A. 3-2) 

©  (A. 3-3) 

(A.3.M 

(A.3.5) 

(A. 3*6) 

(A. 3*7) 


Integrating  (A. 3- 7): 
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(A. 3.8) 


Substituting  for  £e  from  (A. 3*1) 


£  -4- 


«o  60  \  f 

o©+?Zvv  >  p6*-* dr  (A,39) 

K  w 


k  >  I 


MS  I 


Conparing  coefficients  of  K, 


(A. 3. 10) 


'V 


(A. 3- 11) 


Also  from  (A.3.1) > 


_  G(-,tO  *(^) 


To  determine 


5£© 


we  make  use  of  the  expression 


(A. 3.12) 


O 
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_  1  S'  4  -  1&&  \  (a.3.13) 

ir4  ~  ^  L  J 


which  in  view  of  (A.3-U)  “d  (A-3-6)  beccanes 


>te  ,  ^  ©+  -  Y  -  &  *0  ^ 

—  ■  0  *  ^  ur- 

h=.  1  r\ 


re*-'df 

0-.1  "v 


In  view  of  (A.3.1)  and  equating  powers  of  K: 


(A.3.1M 


,*1.\  =  (#-3'15) 


(*•)  = 
v-  >*■ 

^  r  r  'v.  dr 

(A.  3-16) 

Hence  from  (A. 3*12) 

<2 

(A. 3- 17) 

•Y  lo-  — 

By  definition  and  In 

view  of  (A.3.1), 

(A.3.18) 


l6l 


where 


S_  r  Gfr,  0  *t,KKO 


(A.  3- 19) 


It  can  be  readily  shown  that 


2  £  \ 

3  1  9«v  •  J 


Therefore  in  view  of  (A. 3*19)  and  (A. 3*20), 


(A.3-20) 


S  -  i  &  *  ?  2^(|i8)  C  l  (A. 3.21) 


and  from  (A. 3-1 8) 


Ir&VJ 


(A.3.22) 


Eq.  (A. 3*22)  formally  completes  the  cycle  of  operations  in  as  far 
as,  is  found  from  (A.3»10),  (?£.•)  from  (A.3.15),  and  £-r 

o  •j'f  /©  0 

from  (A. 3. 17).  Hence  <o0  is  found  from  (A. 3-22),  and  consequently 

£©  and  (^L®)  can  be  found  from  (A. 3 *11)  and  (A. 3*16)  respec- 
'  ' 1 

tlvely.  This  procedure  can  be  repeated  for  higher  values  of  M,  . 


Criteria  of  Convergence 


From  (A. 3. 11),  (A. 3-16)  and  (A. 3.22)  we  obtain 
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u-3-23) 


also  from  (A. 3-23) 


i<*js 


A 

\G*  Hk5,"'  ^  j +  [  £*{"  hk.,+ i  ‘v- 


(A. 3-24) 


where 


H 


V*  - 1 


‘fe’-Mr*--  dr 


(A. 3.24a) 


is  a  non- negative  monotonically  increasing 
function  of  ^  for  all  n  and  "t 
Hence* 


K  6ft,  *,*1  dt 

<  &  I  H  ™ 

°  3* 

0  1  a  - 1 

(A.3.25) 


where 


Also 


I  K.-H  -  K-.L*  t- 

« '  .  «  « 

^  &  *  H„_,  ch  J  ^  [  |  6  * H a-.  I cN  <  ^GJhJ^|ck  (A.3.26) 


For  proof  see  end  of  Appendix. 
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in  view  of  (A.3.25). 
However 


Then  froa  (A. 3*28),  (A. 3-25)  (A. 3*2^) 


We  now  write 


* 


H 


Integrating  by  parte 


introducing  (A-3-31)  ^  U-3.29)  we  finally  obtain 


or 


(A. 3-26) 


(A.3.29) 


(A. 3- 30) 


(A.  3. 31) 


(A. 3- 32) 


<r,'6' 


(A. 3-33) 


Hence  a  sufficient  (but  not  necessary)  condition  of  convergence  is 
that 


A  typical  value  of  'ixf^  would  be  2.0. 
In  this  case 


Go 

K 

which  in  terms  of  \7m 


becomes, 


3(l-2  O 


> 


O.S 


(A. 3. 34) 


(A. 3. 35) 


(A.3.36) 


It  is  important,  however,  that  smaller  values  of  do  not 

exclude  convergence  although  their  sufficiency  cannot  be  established. 
Nevertheless  we  expect  that  for  most  viscoelastic  materials  v»  will 
be  much  nearer  l/2,  and  hence  (A. 3*36)  1*  sufficient  for  most  practical 


purposes . 
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To  Shew  that  If  Monotonically  IacrtMlm  In  ^ 

-t 

|  ^  GC*  ci^i  |  £  (H  1  Gq  (A. 3-37) 


where 


H  =  H 


wax 


IK 


C  G  (■*,  t,t)  =  G  ( M  -  H  ")  +  G  (H  -  H  )•+ — *  G  (H  - M  ^ 

/  1  ^  |  1  °  '  Z  1  M-H. 

&,  Wo  -*  H,(GI-C1|U  ***  *-^6,  | 

<  <S,|HJ  -(G.-S.^H  |  *  ...  *  GJHJ  ?tIHUSk  <*-3-38) 


where  Q  is  acne  value  of  Q 

on 

In  the  limit  £(„  ^ 


in  the  nth  interval 
GC-v  t  -f)  =  G  >  and  therefore 


<  I H  I  So 


(A.3-39) 
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GENERAL  DISCUSSION 

When  the  critical  dependence  of  the  material  properties  of  visco¬ 
elastic  solids  on  temperature  was  established,  it  was  realized  that  the 
prospects  of  obtaining  analytic  solutions  to  the  related  boundary  value 
problem  with  arbitrary  geometry,  were  rather  poor.  In  fact,  up  to  the 
time  the  present  work  was  undertaken,  the  only  two  problems  that 
yielded  exact  solutions  in  the  formal  sense,  were  those  of  the  infinite 
state  and  the  sphere  with  polar  symmetry,  these  solutions  being  limited 
to  solids  with  thermorheologically  simple  behavior. 

The  purpose  of  the  present  work  was  dual  in  tie  sense  that  though 
the  infinite  cylinder  was  initially  the  central  problem  to  which  a 
solution  was  sought,  general  approximate  techniques  that  would  apply 
to  any  desired  geometry  were  developed,  and  tested  with  very  encourag¬ 
ing  results  in  the  two  cases  where  exact  solutions  were  known. 

The  cylinder  problem  was  also  treated  successfully,  and  two 
analytic  solutions  were  given.  The  perturbution  solution  can  be  viewed 
as  an  "exact"  solution,  in  -s  far  as  the  series  expansion  of  the  unknown 
function  has  been  proved  convergent,  and  hence  the  desired  degree  of 
accuracy  may  be  achieved  by  calculating  a  sufficient  number  of  terms. 

Also,  within  the  limitation  of  material  incompressibility  the 
problems  of  the  sphere  and  cylinder  were  examined  and  solved  in  the 
presence  of  inertia  forces.  It  is  worth  noting  that  in  the  presence  of 
inertia  forces,  compressibility  precludes  closed  form  solutions  even  in 
the  case  of  the  infinite  slab,  since  the  equilibrium  equation  can  no 
longer  be  integrated  directly. 
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A  perturbation  technique  is  indicated  here  aa  veil,  and  will  form 
the  baaia  of  future  vork. 

An  interesting  consequence  of  the  Material  dependence  on  temperature, 
from  a  structural  integrity  design  standpoint,  ia  the  effect  of  a 
periodic  temperature  history  on  the  stresses  in  a  viscoelastic  slab. 

It  was  shown  in  Appendix  I  that  the  maximum  stress  was  almost 
twice  the  value  that  would  have  occurred  in  an  elastic  slab  under  the 
same  temperature  history.  This  is  contrary  to  views  expressed  in  the 
literature  Q13]  ,  that  stresses  due  to  a  step  input  are  likely  to  be  the 
design  stresses. 

In  fact  this  phenomenon  is  only  but  one  facet  of  the  '"Ifcermal 
Cycling  Problem"  which  is  of  particular  concern  in  systems  with  solid 
propellant  configurations.  In  view  of  the  above  result  it  is  natural 
to  inquire  into  the  existence  of  a  critical  thermal  cycle  which  will 
produce  the  maximum  possible  stress  at  some  point  in  the  viscoelastic 
body. 

Though  the  existence  of  such  a  cycle  is  instinctively  certain,  its 
determination  is  a  formidable  problem,  and  it  requires  further  studies. 

Finally  in  the  third  Chapter,  the  horizontal  slump  problem  of  a 
viscoelastic  cylinder  under  isothermal  conditions,  contained  in  an 
elastic  shell  is  solved  formally.  Numerical  computations  are  under 
way. 

This  is  a  problem  of  concern  when  solid  propellant  systems  have 
to  be  stored  for  a  prolonged  period  of  time. 
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